AEROSPACE ENGINEERING 
AE449 Senior Design Project II 
Auburn University, Alabama 


FINAL STUDY REPORT FOR THE 
ADVANCED SPACeFtRANSPORTATION SYSTEM 

Volume II: Technical Report 


A-C r :-lo^ 953 ) i°Ar 
i. r°ANS p QPT AfluN 

rr (Auburn Univ.) 


c jhut ( l i.. ^ « n vbsr: n 

Y " T " M, VUl.UML t i n -' 1 
: L°<CL "V 


- i 


/1 


N U 0 — < *-■ o ‘ ■ ■> 


unci < i s 
0?b4^j 


Submitted to: Dr. James O. Nichols 


Submitted by: 


James N. Adinaro 
Philip A. Benefield 
Shelby D. Johnson 
Lisa K. Knight 


Date Submitted: April 27,1989 



Table of Contents 


1.0 Configuration Determination 

2.0 STME Performance Calculation Summary 

3.0 Mission Analysis 

4.0 Orbital Analysis 

5.0 Stability and Control Analysis 


1 

7 

14 

49 

53 







1.0 Configuration Determination 


To determine the best configuration from all candidate configurations !! was pessary 
first to calculate minimum system weights and performance. To optimize the desgn 
neceaary to vary configuration-specific variables such as total system weight, thrust-to- 
weightTatiosJoum durations, fetal thrust available, and mass fraction for_ the'system* 
Ootimizine each of these variables at the same time is technically unfeasible and not 
necessarily mathematically possible. However, discrete sets rf datat can begenerai e 
which will eliminate many candidate configurations. From the p S 

remaining designs, a final configuration can be selected. 

To ease this optimization process an interactive spreadsheet wascreatedltacceptedinitml 
thrust to weight ratios, number of main engines, vehicle diy weight, a«d nitml (p e- 
seDaration) bum time. From this information, it calculated maximum weights and 
propellant mass fractions available for each vehicle and the final post-separation u 

time for the orbiter. 

After selecting the number of engines and thrust to weight ratios for both thebooster*he 
orbiter the length of the initial bum was varied to maximize fuel available to the orbiter for 
the second phafebum. This maximum weight was then used to determine tank sizings 
which dictated total volume necessary within the booster and orbiter. 

Included are the three most important designs considered: one which closely approximates 
the design criteria set forth in a Marshall Space Flight Center study of the Shuttle II, the 
configuration used in the initial proposal; and the final configuration presented and 
analyzed in this report. Note that these three systems represent important iterative steps in 

our final design. 

A listing by cell of the formulas used to generate the aforementioned data is included for 
reference. 
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ORBITER DATA - 


.Thrust/Weight 

1.2000 

Initial Boost Engines 

4.00 

Total Thrust at Takeoff 

1, 454,100.00 lbf 

Available Takeoff Weight 

1,211,750.00 Ibm 

Dry Weight 

180,000.00 Ibm 

Payload 

60,000.00 Ibm 

RCS Fuel (MMH/N04) 

20,000.00 Ibm 

Total Orbiter Weight 

260, 000.00 Ibm 

Maximum Fuel Weight 

951,75 0.00 Ibm 

Fuel Available Initially 

951,750.00 Ibm 

Initial Boost Engines 

4.00 

Initial Fuel Consumed 

505,080.80 Ibm 

Primary Fuel Available 

448,€89.20 Ibm 

Auxilliary Fuel 

4 9 6, 95 2.57 Ibm 

Total Fuel Available 

94 3, 621.. 7 7 Ibm 

Final Boost Engines 

4.00 

Approximate Final Thrust 

1,740,000.00 lbf 

Approximate Final T/W 

1.4456 

BOOSTER DATA - 


Thrust/Weight 

1.4000 

Initial Boost Engines 

5.00 

Total Thrust at Takeoff 

1, 817,825.00 lbf 

Available Takeoff Weight 

1,298,303.57 Ibm 

Dry Weight 

160,000.00 Ibm 

RCS Fuel (MMH/N04) 

10,000.00 Ibm 

Total Booster Weight 

170,000.00 Ibm 

'Fuel Weight 

1,128,303.57 Ibm 

Fuel Available Initially 

1,128,303.57 Ibm 

Initial Boost Engines 

5.00 

Initial Boost Fuel 

631,351.00 Ibm 

Reservoir Fuel Pumped 

4 96, 95 2.57 Ibm 

OVERALL DATA ------- 


Thrust/Weight 

1.3034 

Initial Boost Engines 

9.00 

Total Thrust at Takeoff 

3,271,725.00 lbf 

Available Takeoff Weight 

2,510,053.57 Ibm 

Dry Weight 

340, 000.00 Ibm 

Payload 

60,000.00 Ibm 

RCS Fuel (MMH/N04) 

30,000.00 Ibm 

Total Dry Weight 

430,000.00 Ibm 

Fuel Weight 

2,080,053.57 Ibm 

Fuel Available Initially 

2,080,053.57 Ibm 

Initial Boost Engines 

9.00 

Initial Boost Fuel 

1,136,431.80 Ibm 

Fuel Available Finally 

943,621.77 Ibm 


Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


133.00 sec 
248.48 sec 
3 81.48 sec 

363 , 525.00 lbf 
435 , 000.00 lbf 

949.40 Ibm/sec 
0.0700 ft~3 
1.1490 ft A 3 
62.43 lbm/ft~3 

135,964.29 Ibm 
815,785.71 Ibm 
31,112 . 40 ft A 3 
11,372.68 ft*3 
42,485.07 ft-3 
88.0301 ft 
21.4203 ft 
109.4504 ft 


161,186.22 Ibm 
967,117.35 Ibm 
36,883.88 ft~3 
13,482.36 ft~3 
50,366.23 ft~3 
69.4704 ft 
25.3939 ft 
94.8643 ft 


V ■ 
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ORBITER DATA - 


Thrust/Weight 

1.2000 

Initial Boost Engines 

5.00 

Total Thrust at Takeoff 

1, 817,825.00 lbf 

Available Takeoff Weight 

1,514,687.50 lbm 

Dry Weight 

180 , 000.00 lbm 

Payload 

60 , 000.00 lbm 

RCS Fuel (MMH/N04) 

20,000.00 lbm 

Total Orbiter Weight 

260 , 000.00 lbm 

Maximum Fuel Weight 

1,254,687.50 lbm 

Fuel Available Initially 

1,254,687.50 lbm 

Initial Boost Engines 

5.00 

Initial Fuel Consumed 

564,893.00 lbm 

Primary Fuel Available 

68 9, 794.50 lbm 

Auxilliary Fuel 

563,410.57 lbm 

Total Fuel Available 

1,253,205.07 lbm 

Final Boost Engines 

4.00 

Approximate Final Thrust 

1,740,000.00 lbf 

Approximate Final T/W 

1.1499 

BOOSTER DATA - 


Thrust/Weight 

1.4000 

Initial Boost Engines 

(_T 

m 

o 

o 

Total Thrust at Takeoff 

1,817,625.00 lbf 

Available Takeoff Weight 

1,298,303.57 lbm 

Dry Weight 

160,000.00 lbm 

RCS Fuel (MMH/N04) 

10,000.00 lbm 

Total Booster Weight 

170,000.00 lbm 

Fuel Weight 

1 , 128 ,303.57 lbm 

Fuel Available Initially 

1,128,303.57 lbm 

Initial Boost Engines 

5.00 

Initial Boost Fuel 

564,893.00 lbm 

Reservoir Fuel Pumped 

563,410.57 lbm 

OVERALL DATA - 


Thrust/Weight 

1.2923 

Initial Boost Engines 

10.00 

Total Thrust at Takeoff 

3,635,250.00 lbf 

Available Takeoff Weight 

2 , 812,991.07 lbm 

Dry Weight 

340,000.00 lbm 

Payload 

60,000.00 lbm 

RCS Fuel (MMH/N04) 

30,000.00 lbm 

Total Dry Weight 

430 , 000.00 lbm 

Fuel Weight 

2,382,991.07 lbm 

Fuel Available Initially 

2,382,991.07 lbm 

Initial Boost Engines 

10.00 

Initial Boost Fuel 

1,129,786.00 lbm 

Fuel Available Finally 

1,253,205.07 lbm 


Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


119.00 sec 
330.00 sec 
44 9.00 sec 

363 , 525.00 lbf 
435 , 000.00 lbf 

949.40 Ibm/sec 
0.0700 ft ^3 
1.1490 ft *3 
62.43 lbm/ft~3 

179,241.07 lbm 
1 ,075,446.43 lbm 
41,015.32 ft A 3 
14 , 992.55 ft * 3 
56,007.87 ft*3 
116.0496 ft 
28.2383 ft 
144.2879 ft 


161,186.22 lbm 
967,117.35 lbm 
36,883.88 ft~3 
13,482.36 ft~3 
50,366.23 ft^3 
69.4704 ft 
25.3939 ft 
94.8643 ft 
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ORBITER DATA - 


Thrust/Weight 

1.2000 

Initial Boost Engines 

4.00 

Total Thrust at Takeoff 

1, 454,100.00 lbf 

Available Takeoff Weight 

1,211,750.00 Ibm 

Dry Weight 

180,000.00 Ibm 

Payload 

60,000.00 Ibm 

RCS Fuel (MMH/N04) 

20,000.00 Ibm 

Total Orbiter Weight 

260,000.00 Ibm 

Maximum Fuel Weight 

951,750.00 Ibm 

Fuel Available Initially 

951,750.00 Ibm 

Initial Boost Engines 

4.00 

Initial Fuel Consumed 

562,044.80 Ibm 

Primary Fuel Available 

389,705.20 Ibm 

Auxilliary Fuel 

544,897.09 Ibm 

Total Fuel Available 

934, 602.29 Ibm 

Final Boost Engines 

4.00 

Approximate Final Thrust 

1, 740,000.00 lbf 

Approximate Final T/W 

1.4544 

BOOSTER DATA ------- 


Thrust/Weight 

1.4000 

Initial Boost Engines 

6.00 

Total Thrust at Takeoff 

2,181,150.00 lbf 

Available Takeoff Weight 

1,557,944.29 Ibm 

Dry Weight 

140,000.00 Ibm 

RCS Fuel (MMH/N04) 

10,000.00 Ibm 

Total Booster Weight 

170,000.00 Ibm 

Fuel Weight 

1, 387,944.29 Ibm 

Fuel Available Initially 

1,387,944.29 Ibm 

Initial Boost Engines 

6.00 

Initial Boost Fuel 

843,047.20 Ibm 

Reservoir Fuel Pumped 

544,897.09 Ibm 

OVERALL DATA ------- 


Thrust/Weight 

1.3125 

Initial Boost Engines 

10.00 

Total Thrust at Takeoff 

3,435,250.00 lbf 

Available Takeoff Weight 

2,749,714.29 Ibm 

Dry Weight 

340,000.00 Ibm 

Payload 

40,000.00 Ibm 

RCS Fuel (MMH/N04) 

30,000.00 Ibm 

Total Dry Weight 

430,000.00 Ibm 

Fuel Weight 

2,339,714.29 Ibm 

Fuel Available Initially 

2,339,714.29 Ibm 

Initial Boost Engines 

10.00 

Initial Boost Fuel 

1, 405,112.00 Ibm 

Fuel Available Finally 

934,402.29 Ibm 


Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


148.00 sec 

246.10 sec 

394.10 sec 

363 , 525.00 lbf 
435 , 000.00 lbf 

949.40 lbm/sec 
0.0700 ft *3 
1.1490 ft ^3 
62.43 lbm/ft-3 

135,964.29 Ibm 
815,785.71 ibm 
31,112.40 ft^3 
11,372.68 ft*3 
42,485.07 ft~3 
88.0301 ft 
21.4203 ft 
109.4504 ft 


198,280.61 Ibm 
1,189,683.67 Ibm 

45.372.10 ft*3 

16.585.10 ft~3 
61, 957.20 ft^3 

85.4579 ft 
31.237 9 ft 
116.6958 ft 
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^ > ORBITER DATA - 

Thrust/Weight 
Initial Boost Engines 
Total Thrust at Takeoff 
Available Takeoff Weight 
Dry Weight 
Payload 

RCS Fuel (MMH/N04) 

Total Orbiter Weight 
Maximum Fuel Weight 

Fuel Available Initially 
Initial Boost Engines 
Initial Fuel Consumed 
Primary Fuel Available 
Auxilliary Fuel 
Total Fuel Available 
Final Boost Engines 
Approximate Final Thrust 
Approximate Final T/W 

BOOSTER DATA - 

Thrust/Weight 
Initial Boost Engines 
Total Thrust at Takeoff 
Available Takeoff Weight 
Dry Weight 
RCS Fuel (MMH/N04) 

Total Booster Weight 
^ Fuel Weight 

Fuel Available Initially 
Initial Boost Engines 
Initial Boost Fuel 
Reservoir Fuel Pumped 

OVERALL DATA - 

Thrust/Weight 
Initial Boost Engines 
Total Thrust at Takeoff 
Available Takeoff Weight 
Dry Weight 
Payload 

RCS Fuel (MMH/N04) 

Total Dry Weight 
Fuel Weight 


1.2 

5 

-B3*E6 

-B4/B2 

180000 

60000 

20000 

-B6+B7+B8 

-B5-B9 

-B10 

-B3 

-B13 *18*X2 
-B12-B14 

— B3 4 

—B15+B16 
4 

-B18*X7 

— B19/(B9+B17) 


1.4 

5 

—B2 3 * E 6 
-B24/B22 

160000 

10000 

-B26+B27 

-B25-B28 

— B2 9 

-B23 

— B3 2 *E8 *X2 
-B31-B33 

— ( B4 +B24)/ (B5+B25) 
-B3+B23 

-B4+B24 

-B38/B36 

-B6+B26 

— B7 

-B8+B27 
—B40+B41+B42 
-B10+B29 


Ibf 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 

lbf 


lbf 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 


lbf 

lbm 

lbm 

lbm 

lbm 

lbm 

lbm 


Fuel Available Initially 
Initial Boost Engines 
Initial Boost Fuel 
Fuel Available Finally 


-B12 +B31 

— B37 

— B4 7 *E2*E8 
-B46-B48 


lbm 

lbm 

lbm 
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Initial Boost Burn Time 
Final Boost Burn Time 
Total Burn Time 

Tsl 

Tvac 

Wdot 

SG - Hydrogen 
SG - Oxygen 
Density of Water 

Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


119 

-B17/(318*18) 
-X2+X3 

363525 

435000 

949.4 

0.07 

1.149 

62.43 

-B12/7 

-B12*6/7 

«X13/(X9*X11) 

-X14/ (X10*X11) 
-X15+X16 

-X15/(PI()/4*15 A 2) 

• 116/(PI()/4*2 6 A 2 ) 

• XI0 +X19 


sec 

sec 

sec 

lbf 

Ibf 

lbm/sec 
ft A 3 
ft A 3 

lbm/ft*3 

Ibm 

lbm 

ft~3 

ft A 3 

ft~3 

/ 2 ft 

ft 
ft 


Max Fuel - Hydrogen 
Max Fuel - Oxygen 
Volume - Hydrogen 
Volume - Oxygen 
Total Volume 
Eq. Hydrogen Tanks 
Eq. Oxygen Tanks 
Minimum Length Necessary 


B29/7 

lbm 

B29*C/7 

lbm 

127/(19**11) 

ft *3 

128/(110*111) 

ft A 3 

129+130 

ft ~3 

129/(?Z()/4 *2 8^2) 

ft 

130/(PX ()/ 4*2 6 A 2) 

ft 

132+133 

ft 
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2.0 STME Performance Calculation Summary 

A model of the total thrust available from the Space Transportation Main Engine 3 used on 
the Shuttle II was necessary to evaluate its launch profile. Thrust is a function of y 
f^abtes mclS^Sexpansion ratio, exit velocity, ambient pressure and exit 
oressure Since the STME uses a translatable skirt, the expansion ratio of the engine 

be set at discrete values of35:1,70:1, and 150:1. To insure that th ® e "^ 
efficiently during all points in the launch, it was necessaiy to calculate the altitudes at 
which the skirt should be translated. This analysis was performed by calculating thrust 
coefficients for all three nozzles and finding the points of intersection. 

These calculations were based on the data presented in the Aeroject contractor's report er 
ffie 4e ?962 standard atmosphere was used as a model ^pressured density as 

functions of altitude. Most of the engine data was solved using simple, iterative 
techniques, isentropic relations, and graphical methods. 
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2.1 Analysis 

Given an expansion ratio for a particular nozzle and an appropriate value for the 
coefficient of specific heats, the exit Mach number may be determined. 

Assuming y = 1.24 for a mixture of liquid oxygen and liquid hydrogen in a 6:1 mass ratio 
and the first (fixed) nozzle expansion ratio of 35:1, this isentropic relation holds: 

Ac - J— | i i77*»i 

/\* *1«-L J 

This function may be solved iteratively to obtain exit Mach numbers for the nozzles. This 
was solved using the included iterative, and the following were obtained as the exit Mach 

numbers: 




l« Expansion Ratio 
35:1 
70:1 
150:1 


Frit Maril Number 
4.369 
4.931 
5.581 


Next, the flow in a nozzle may be assumed to be very nearly isentropic, such that the 
chamber pressure may be related to the nozzle exit pressure through the isentropic relation 

^=r - ( 1 ♦ *zf 

Given a c hamb er pressure of 2385 psia under normal power levels and 3125 psia under 
emergency conditions, the following exit pressures were obtained. 


Pe (NFL) PefEPL) 
For e = 35:1 5.069 6.642 

For e = 70:1 2.058 2.696 

For e = 150:1 0.771 1.010 


For a given expansion ratio, a nozzle will have a maximum efficiency at a certain design 
altitude. At this altitude, the exit pressure is equal to the ambient pressure of the 
atmosphere. From the 1962 standard atmosphere, the design altitudes were determined to be 

hd (Normal Power Level) (Emergency Power Levels) 

For e = 35:1 26,700 ft 20,400 ft 

For e = 70:1 45,800 ft 40.200 ft 

For e = 150:1 66,300 ft 60,600 ft 


In order to obtain the proper altitudes at which the sections of the nozzle should be extended, 
the following procedure is followed. 


The thrust coefficient is defined as: 

Or = 


Pr A* 


Since thrust varies with altutde, this expression may be rewritten as. 


c 



Z y* 


Y“ I / \ r 




fcj" 1 


Pr 


Pc_ 

Pr 


PA 


Ac 

A * 
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For 


example, consider the nozzle with an expansion ratio of 70:1 at an altitude of 50,000 feet. 


/ vP /«•«SaV ,1aM 1 t 

C T s ‘ fed?) J CS&-I 

Similarly, thrust coefficients for other nozzle and atmospheric conditions maybe 
calculated. As shown in Figure 2.1.1, the points where these graphs cross i indicate 
altitudes where the nozzles should be translated for maximum efficie cy. 

For Normal Power Levels, the translation altitudes were found to ^36,900 feet and 56,500 
feet For Emergency Power Levels, the altitudes were 31,100 feet and 50,800 feet. 


Z.O60 * /*6.2>w 
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11000 
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13000 


14000 


15000 


16000 


17000 


18000 


19000 


20000 


21000 


22000 


23000 


24000 


25000 


26000 


27000 


28000 


29000 


30000 


31000 


32000 


33000 


34000 


35000 


36000 


37000 


38000 


39000 


40000 


41000 


42000 


43000 


44000 


45000 


46000 


47000 


48000 


49000 


50000 


51000 


52000 


53000 


54000 


55000 


14.7000 


14.1700 


13.6600 


13.1700 


12.2300 


11.7800 


10.1100 


7.6470 


7.3390 


5.6960 


5.2200 


4.364 


3.6260 


3.142 


2.7200 


2.4710 


1.9430 


1.8520 


1.7650 


1.6068 


1.7840 


1.8052 


1.8149 
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Figure 2.1.1: Plot of Thust Coefficients vs. Altitude (NPL) 
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Figure 2.1.2: Thust Coefficients vs. Altitude (EPL) 
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3.0 Mission Analysis 

To prove the flight worthiness and technical feasibility of the ASTS as proposed, a launch 
profile was constructed and simulated. This simulation consisited of a dynamic model of 
the vehicle’s instantaneous position, velocity, acceleration, thrust, weight, orientation, 
and flight path with respect to the horizon. In addition to this vehicle-specific data, 
atmospheric conditions such as temperature, pressure, and density were modeled. 


The basic format of this simulation was a fourth-order Runge-Kutta integration scheme. A 
number of asumptions were necessary to develop this model of the launch trajectory. First, 
the model is constrained by a run time of approximate seven minutes and a total distance 
downrange not exceeding 800 nautical miles, so the surface of the earth was modeled as 
being flat. In addition, rotational effects of the earth were neglected. This assumption 
cannot degrade the validity of results, since it could only create conditions more favorable 
to a due east launch from Kennedy Space Center. Drag forces were neglected for two 
principal reasons: they can be considered small in comparison to thrust forces necessary 
for a ballistic trajectory, and they would not exceed the advantages gained by 
incorporating the earth's rotation. 

All initial conditions were set based on values calculated from the design process outlined 
in Section 1.0. Bum times, initial weights, total thrust available, and fuel available are the 
the most noteworthy examples. Thrust was computed as a funciton of altitude, based on 
translation altitudes for the extendable skirt on the STME nozzles, pressure forces, and 
known values of sea level and vacuum values. It was estimated that the maximum 
acceptable pitch rate was 0.5 degrees per second to retain stability during launch. This is 
the value used in the program for all pitch manuevers. 

Before finalizing the launch profile, it was necessary to establish the range of orbits into 
which the orbiter could be successfully launched. Based on current STS mission analyses 
and estimations of maximum altitude and velocity available at this altitude, our launch 
goal was projected to be an altitude of 100 nautical miles and a velocity equal to that of the 
local circular speed. Inclination was varied during flight both to achieve this orbit and to 
place the booster in an alignment practical for separation. 








program AE447Prog; 


($MC68020+} 

{$MC68881+} 

{SSETC Elems881:-true) 


Written by Del Johnson (in TML Pascal II vl-00 under MPW v2. 02)______ 

Written””" AE447-AE449 Shuttle II design g roup (Fall 1988 - s P rir ^_j;f^_ 

;"rp''g"«'irdeII^d‘"“";r"r""hrof d,,i,n configuration^for the 

execution of these programs was used in an iterative process to fine-tune the 
model as presented in the report. _ 

Definition of variables u[i]: 


( Integral variables.- } 


u [ 1] 
u( 2] 
u[ 3] 
u [ 4] 
u[ 5] 
u[ 6] 


u[x] - Location along the x-axis 

u[xd] - Time rate of change of location along the x-axxs 
u[v] * Location along the y-axis 

u[yd] * Time rate of change of location along the y-axis 
ufth] - Angle with respect to the y-axis 

utthd] - Time rate of change of angle with respect to the y-axis 


{ Other variables— } 


u[alpha]- Angle of attack 
u[gamma]- Flight path angle 
u[wt 3 - Total vehicle weight 

u[wtd] - Time rate of change of total vehicle weight 
u[dens] - Density of the atmosphere 
u[press] - Density of the atmosphere 
u[ thrust] - Thrust produced by all engines firing 


uses MemTypes,SANE; 


const 

{ Number of equations- } 

numlntEquations - 6; 

numTotEquations — 11; 

{ Time Considerations- } 

tbegin - 0.00; { seconds } 

tend - 420.00; { seconds } 

dt - 0.10; { seconds } 

cint * 10; 

cols - 4; 

{ Subscripts of the variables of motion— } 
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xd ■ i ; 

x " 2; 

yd ■ 3 ; 

y ■ 4; 

wt “ 5; 

th - 6; 

vtot " 7; 

gamma ” 8; 

thrust “ 9; 

press * 10; 

dens *■ 11; 

Separation “ 1 
Shutdown “ 2 


StartPitchl “ 3 
StopPitchl ■ 4 
StartPitch2 “ 5 
StopPitch2 • 6 


l Constants related to the problem™ } 

{ Mathematical constants™ } 

conv - 57.29577951; { degrees/radians } 


{ 


Physical 

Re 

gO 

rho 

rhoO 

Psl 


constants™ } 

- 20938912.00000; ( feet } 

- 32.17400; { feet/second~2 } 

- 0.0023768; { pound-mass/feet^ } 

- 0.0023768; < pound-mass/feet''3 } 

- 14.70000; l pound-force/inchest } 


{ Vehicle characteristics™ ) 

thrO - 363535.00000; { lbf } 

wtdO - -949.40000; { lbm/sec } 


pressAlt “ 25000 
densAlt ” 25000 

Ael - 3247 
Ae2 - 6503 
Ae3 - 14526 


00000; { feet } 
00000; { feet } 

22085; { inches A 2 } 
88219; { inches A 2 } 
72443; { inchest } 


designAlt1NPL 
designAlt2NPL 
designAlt3NPL 
designAlt1EPL 
designAlt2EPL 
designAlt3EPL 


26700.00000; { feet } 
45800.00000; { feet } 
66300.00000; { feet } 
20400.00000; { feet } 
40200.00000; { feet } 
60600.00000; { feet l 


transAltlNPL - 36900.0000; { feet } 
transAlt2NPL - 56500.0000; { feet } 
transAltlEPL - 31100.0000; { feet } 
transAlt2EPL - 50800.0000; { feet J 


PelNPL - 5.069; { lbf/feet A 2 } 
Pe2NPL - 2.058; { lbf/feet^ } 
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Pe3NPL - 0.771 
PelEPL - 6.642 
Pe2EPL - 6.296 
Pe3EPL - 1.010 

thrMaxlNPL - 4 
thrMax2NPL - 4 
thrMax3NPL “ 4 

thrMaxlEPL « 5 
thrMax2EPL - 5 
thrMax3EPL ■ 5 


1; { lbf/feet A 2 } 

2; { lbf/feet A 2 } 

6; { lbf/feet^2 } 
0; { lbf/feet*2 } 

411259.15000; 

424577.04000; 

435000.00000; 

545789.15000; 

556312.08000; 

580000.00000; 


{ 


Problem intial conditions... } 


xdO 


0.00; 

t 

feet/second 

> 

xO 

- 

0.00; 

{ 

feet ) 


ydO 

- 

0.00; 

{ 

feet/second 


y0 

- 

0.00; 

{ 

feet } 


thO 

- 

90.00; 

{ 

degrees } 


vtotO 

« 

0.00; 

{ 

feet/second 

} 


type 

uArray 
name Array 
printrange 


array 
array 
1..32; 


[1..numTotEquations] 
[1..numTotEquations] 


of extended; 
of str255; 


var 


outfile 
tab 
3p 


text; 

char; 

array 


[1..numTotEquations] of integer; 


DryWtBooster, DryWtOrbiter, wtO : extended; 


extended; 


time 

k 

u,uu 
if j/m r n 
printset 


array [1..6] of extended; 

array [1..4 f 1..numlntEquations] of extended; 

uArray; 

integer; 

set of printrange; 


{ Variables of motion^. } 

wtd, thr, thd : extended; 

thrMax : extended; 


function g(alt : extended): extended; 
begin 

g g0*sqr (Re)/sqr (Re+alt) ; 


end; 

function density(alt : extended): extended; 
begin 

density :* rho0*exp(—alt/densAlt); 


end; 

function pressure(alt : extended): extended; 
begin 

pressure :- Psl*exp(-alt/pressAlt); 

end; 



17 










function pitchRate(t : extended): extended; 
begin 

if (t < Time[StartPitchl]) then pitchRate 0.00 else 

if (t < Time[StopPitchl] ) then pitchRate -0.50 else 

if (t < Time[StartPitch2]) then pitchRate :* 0.00 else 

if (t < Time[StopPitch2] ) then pitchRate -0.50 else 

if (t < tend) then pitchRate 0.00; 

end; 


function thrustAvail(alt : extended): extended; 
var 

Ae,thrl : extended; 

begin 

if (alt<transAltlNPL) then begin 
Ae Ael; 
thrl :» thrMaxlNPL; 

end 

else 

if (alt<transAlt2NPL) then begin 
Ae Ae2; 
thrl :* thrMax2NPL; 

end 

else begin 

Ae :* Ae3; 

thrl :* thrMax3NPL; 

end; 

thrustAvail :** thrl - pressure(alt)*Ae; 

end; 

procedure initVariables; 
begin 

t tbegin; 
u[yd] ydO; 
u [y] yO; 

u[xd] xdO; 
u[x] xO; 

u[wt] :* wtO; 
u(th] :* thO; 
u[gamma] :* 90; 

u[vtot] sqrt(sqr(u[yd])+sqr(u[xd])); 
u[press] :« pressure(u[y]); 
u[dens] :» density(u[y]); 

end; 

procedure print(style : integer); 
var 

i : integer; 
begin 

case style of 
1 : begin 

write(outfile, f t ■ f ,t:6:2); 

for i :■ 1 to numTotEquations do begin 

if (i in printset) then write(outfile, 1 u[’,i:2, f ] 

-',u[i]:12:2); 

if (i in printset) then write( 1 u[',i:2r f ] 

- 1 , u ti] : 12 :2) ; 

end; 
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writeln(outfile); 
writeln; 
end; 

2 : begin 

i :■ 0; 

writeln(outfile,* Time • •,t:7:2); 
writeln ( , Tiine * 1 ,t:7:2); 

for j :■ 1 to numTotEquations do 
if (j in printset) then begin 
i :■ i+1; 

write (out file, 1 u [*, j :2, ' ] ■ ,l , u [ j] : 12 : 2, 1 )/ 

write( 1 ■ B, f u[j]:12:2 r * )> 

if (i mod cols - 0) then begin 
i :« 0; 

writeln(outfile); 
writeln; 
end; 
end; 

if (i<>0) then begin 
writeln(outfile); 
writeln; 
end; 
end; 

3 : begin 

write(outfile,t:7 :2, tab); 

for j :■ 1 to numTotEquations do 

if (j in printSet) then write(outfile,u(j]:12:sp(j],tab); 

writeln(outfile); 
end; 
otherwise; 
end; (case) 
end; 

{ Total function defining variables in terms of one another... ) 
function f(j: integer; t: extended; u: uArray): extended; 

begin 


case j of 
xd 


- thr*cos(u[th]/conv)*g(u[y])/u[wt]; 

» u[xd]; 

- thr*sin(u[th]/conv)*g(u[y])/u(wt] - g(u[y]); 

- u [yd] ; 

- wtd; 

- thd; 


otherwise; 


end; 


end; 


begin 


Textbook(nil) ; 
tab :« chr(09); 

for j :■ 1 to numTotEquations do 
sp[j3 :« 2; 
sp[press] :* 8; 
sp[dens] :* 8; 


open (outfile, 1 l/AE448prog.dat ); 
rewrite(outfile); 
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m numlntEquations; { Number of integration equations } 

n trunc((tend-tbegin)/dt); { Number of steps } 

{ Set names of all variables- } 

printset [xd,x,yd,y,wt,vtot,th,gamma,thrust,press,dens]; 

{ Set design variables— } 

wtO 2769714.29; { lbf } 

DryWtBooster :“ 170000.00; { lbf } 

DryWtOrbiter : m 260000.00; { lbf } 

Time[Separation] :■ 148.00; { seconds } 

Time[Shutdown] :■ 394.00; { seconds } 

Time[StartPitchl] 20.00; { seconds } 

Time[StopPitchl] 140.00; { seconds ) 

Time[StartPitch2] 160.00; [ seconds } 

Time[StopPitch2] Time[StartPitch2] +38.00; [ seconds } 

{ Set initial conditions of all variables- ) 
initVariables; 

{ Echo initial conditions of all variables in the printing set— < 
print(3); 

{ Beginning of time-step loop— } 
for i 1 to n do begin 

{ Computation of other* non time—dependent variables— } 

if (t — Time [Separation]) then u[wt] :“ u[wt] - DryWtBooster; 

thrMax :« thrustAvail(u[y]); 

if (t < Time[Separation]) then thr :- 10*thrMax else 
if (t < Time[Shutdown] ) then thr 4*thrMax else 

if (t < tend) then thr 0*thrMax; 

u[press] pressure(u[y]); 

u(dens] density(u[y]); 

u[thrust] :■ thr; 

if (t < Time [Separation]) then wtd :■ 10*wtd0 else 

if (t < Time[Shutdown] ) then wtd 4*wtd0 else 

if (t < tend) then wtd :« 0; 

thd pitchRate(t); 

if <u[th] > 360.00) then u[th] u[th] - 360.00; 

if (u[thj < -360.00) then u[th] u[th] + 360.00; 

( Runge-Kutta 4th-Order integration loop of all variables- } 

for j :■ 1 to « do 

k[l,j] dt*f(j,t,u); 

for j 1 to m do 

uu(j] :■ u(j] + k(l,j]/2; 

for j :« 1 to m do 

k [2, j] dt*f(j,t+dt/2,uu); 
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u[gamma] 


for j 
for j 
for j 
for j 
for j 


:■ 1 to a do 

uu[j] :” u[j] + k[2,j]/2; 

: - 1 to a do 

]c[3, j] dt*f(j,t+dt/2,uu); 

:« 1 to m do 

uu[ j] :“ u[j] + It 13 , j]; 

1 to a do 

k[4,j] dt*f(j,t+dt,uu); 

:■ 1 to m do .. 

U[ j] u[j] + (k[l, j]+2*Jc[2, j]+2*k[3, j]+k[4,D])/6; 


t :■ tbegin + i*dt; 

u[vtot] :■ sqrt(sqr(u[xd])+sqr(u[yd])); 

if (u[xd] > 0) then u[gamma] arctan (u[yd]/u[xd])*conv 

90.00; 


{ Printing parameters.- } 

if (i mod cint - 0) and (printset <> []) then print(3); 


end; 

close(outfile); 


end. 


else 
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'll 


•in 



210.00 



541043.75 


548643.88 



564007.36 

7791.21 

571771.13 

7846.271 

KEEHTrh 

7901.53! 

587463.73 


348507.33 


350691.16 


2091.68 


2071.32 


1156626.29 


1152828.69 


1149031.09 


1145233.49 


1141435.89 


1137638.29 


1133840.69 


1130043.09 


1126245.49 


1122447.89 


1118650.29 


1114852.69 


1111055.09 


1107257.49 


1103459.89 


1099662.29 


1095864.69 


1092067.09 


1088269.49 


1084471.89 


1080674.29 


1076876.69 


1073079.09 


1069281.49 


1065483.89 


1061686.29 


1057888.69 


1054091.09 


1050293.49 


1046495.89 


1042698.29 


1038900.69 


1035103.09 


1031305.49 


1027507.89 


1023710.29 


1019912.69 


1016115.09 


1012317.49 


1008519.89 


1004722.29 


1000924.69 


997127.09 


993329.49 


989531.89 


985734.29 


981936.69 


978139.09 


974341.49 


970543.89 


966746.29 


962948.69 


959151.09 































































































































































































l/AE448data 



218 


FT] 

CJI 

m 


ai 



223 


225 


226 


228 


229 


23 


23 


232 


233 


235 


236 


238 


239 


242 


248 


25 


25 


252 


253 


255 


256 


258 


259 


26 


26 


262 


265 


227.00 


228.00 


HI 


230.00 


231.00 


HI 


233.00 


234.00 


235.00 


236.00 


HI 


238.00 


239.00 


240.00 


241.00 


242.00 


243.00 


244.00 


245.00 


246.00 


247.00 


248.00 


249.00 


250.00 


251.00 


252.00 


254.00 


255.00 


256.00 


257.00 


258.00 


259.00 


260.00 


261.00 


262.00 


263.00 


■Kivinn 

8124.69 


8181.01 

216.00 



8294.31 

218.00 


219.00 

8408.50 

220.00 

8465.93 


8523.59 


8581.47 


8639.59 

224.00 

8697.94 

225.00 

■EEmBil 


8815.35 


8874.41 


8933.72 


603377.80 


611418.41 


627667.86 


635877.12 


644143.03 


652465.81 


660845.69 


669282.88 


677777.62 


686330.13 


694940.64 


703609.39 


712336.60 


:lik%VHri 


9113.10 


9233.94 


729967.38 


738871.42 


747834.89 


756858.03 


765941.09 


775084.32 


784287.97 


9417.12 


802877.54 


9540.54 


9602.66 


9665.04 


9727.69 


9790.62 


9853.83 


9917.31 


821711.86 




840793.04 


850426.86 


860123.21 


• I : 


10045.14 


10109.48 


10174.11 


10239.04 


10304.27 


10369.79 


10435.62 


10501.75 


10568.19 


10634.94 


10702.01 


10769.39 


10837.10 


10905.13 


10973.49 


11042.18 


11111.20 


2051.02 


2030.76 


2010.54 


1950.19 


1930.17 


1890.28 


1870.40 


1850.58 


1830.81 


1811.08 


1791.41 


1771.79 


365404.90 


369426.01 


371406.32 



379127.15 


381007.49 


382867.97 


384708.66 


1732.70 


1713.23 


1693.81 


1674.45 


388330.85 


390112.44 


391874.44 


1616.68 


1578.44 


1559.41 


1540.42 


1521.50 


1502.63 


1483.81 


1465.06 


1446.36 


397043.36 


398727.49 


400392.28 


402037.79 


403664.07 


405271.18 


406859.16 


408428.08 


409977.99 


411508.95 


413021.01 


414514.22 


415988.65 


417444.36 


955353.49 


951555.89 


947758.29 


943960.69 


940163.09 


936365.49 


932567.89 


928770.29 


924972.69 


921175.09 


917377.49 


913579.89 


909782.29 


905984.69 


902187.09 


898389.49 


894591.89 


890794.29 


886996.69 


883199.09 


879401.49 


875603.89 


871806.29 


868008.69 


864211.09 


860413.49 


856615.89 


852818.29 


849020.69 


845223.09 


841425.49 


837627.89 


879704.54 

1427.72 

418881.39 

833830.29 

889590.09 

1409.14 

420299.81 

830032.69 

899539.26 

1390.61 

421699.68 

826235.09 


1372.15 

423081.06 

822437.49 

919629.64 

1353.75 

424444.00 

818639.89 

929771.41 

1335.40 

425788.57 

814842.29 

939977.96 

1317.12 

427114.83 

811044.69 

950249.59 

1298.90 

428422.84 

807247.09 

960586.59 

1280.74 

429712.66 

803449.49 

970989.27 

1262.65 

430984.35 

799651.89 

981457.93 

1244.62 

432237.98 

795854.29 

991992.87 

1226 . 65 ! 

433473.60 

792056.69 

1002594.41 

1208.74 

434691.29 

788259.09 

1013262.86 

1190.90 


784461.49 

1023998.53 


437073.12 

780663.89 

1034801.75 

1155.42 

438237.38 

776866.29 

1045672.84 

1137.77 

439383.97 

773068.69 

1056612.12 

1120.20 

440512.95 

769271.09 

1067619.92 

1102.69 

441624.39 

765473.49 

1078696.58 

1085.25 

442718.35 

761675.89 


1067.88 

443794.91 

757878.29 
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l / AE 448 data 



11250.27 


11320.31 


11390.71 


11461.46 


1101057.82 


1123698.561 1016.18 



11604.03 


11748.06 


11820.62 


11893.57 


11966.89 


12040.60 


1146621.60 


1158189.87 


1169829.78 


1181541.71 




1205183.08 


1217113.28 


1229116.99 


1241194.60 



1277874.76 


1290251.92 


12491.141 130270 4.98 



12644.591 1327840.43 


1340523.66 


12877.961 1366123.29 


1379040.54 


1392036.68 


1405112.15 


1418267.40 


1431502.90 


1444819.09 


13438.111 1458216.46 


1471695.48 


13602.37 


1498900.37 


13768.57 


14021.66 


14107.05 


14192.97 


14279.42 1596761.86 


14366 . 4211611084.74 


1625494.88 



14630.70 


14719.93 


14809.72 1684019.21 


14900.11 1698874.07 


14991.08 1713819.62 


15082.651 1728856.44 


15174.841 1743985.13 


1759206.31 


1774520.60 


1789928.63 


444854.13 


445896.09 


446920.84 


447928.48 


448919.06 


449892.65 


450849.34 


451789.20 


453618.71 


454508.51 






457079.04 


457903.18 


458711.11 




460278.64 


461038.41 


461782.28 


462510.36 

_ 


463919.43 


464600.60 



:!•] 


467171.59 


467776.37 


468366.14 


468941.00 


469501.04 


470046.36 


470577.06 


471093.24 


471594.99 


472082.41 


• II 




(Hi: 



:!•] 


473014.70 


473459.76 


473890.92 


474308.27 


474711.93 


475102.01 


475478.61 


475841.85 


476191.83 


476528.69 


476852.53 


477163.47 


477461.64 


477747.14 


478020.11 


754080.69 


750283.09 


746485.49 


742687.89 


738890.29 


735092.69 


731295.09 


727497.49 


723699.89 


719902.29 


716104.69 


712307.09 


708509.49 


704711.89 


700914.29 


697116.69 


693319.09 


689521.49 


685723.89 


681926.29 


678128.69 


674331.09 


670533.49 


666735.89 


662938.29 


659140.69 


655343.09 


651545.49 


647747.89 


643950.29 


640152.69 


636355.09 


632557.49 


628759.89 


624962.29 


621164.69 


617367.09 


613569.49 


609771.89 


605974.29 


602176.69 


598379.09 


594581.49 


590783.89 


586986.29 


583188.69 


579391.09 


575593.49 


571795.89 


567998.29 


564200.69 


560403.09 


556605.49 








































































































































































l/AE448data 




32 


322 


323 


32 


325 


326 


328 


329 


33 


33 


332 


333 


335 


336 


338 


339 


15549.801 1805431^03 


1821028.45 


15741.151 1836721.54 




nil iMiii i f n 

1868397.42 


16033.22 

1884381.57 



1900464.10 


16231.41 

1916645.72 


16331.57 

E&EeEEKEI 


■hwtb 

1949309.11 


RTXSEEMc!c]l 

328.00 

16636.49 

1982377.57 


16739.64 

1999065.57 


16843.56 

2015857.11 

331.00 

16948.27 

iM.M'kUftld 


17053.78 

2049753.92 


17160.10 

2066860.79 


17267.24 

2084074.38 



• H 


336.001 17484.04 



17928.11 


18041.39 


18155.61 


18270.78 


18386.91 


18504.03 




2136363.90 


2154012.84 


2171772.78 


2189644.63 


2207629.31 


2225727.73 


2243940.84 


ME 


2280714.99 


18741.29 


18861.47 


18982.70 


19105.00 


19228.39 


19352.90 


19478.54 


19605.33 


19733.29 


19862.46 


19992.84 


20124.46 


20257.35 


20391.53 


20527.03 


20663.87 


20802.07 


20941.68 


21082.71 


21225.20 


21514.65 


21661.69 


2317959.63 


2336760.93 


* 


2393893.28 


2413183.83 


2432599.45 


2452141.29 


2471810.50 


2491608.27 


2572109.49 


2613163.99 




2654768.60 


2675780.68 


2696934.51 


2718231.56 




•m 



2761261.39 


478280.68 


478528.96 


478765.08 


478989.18 


479201.39 


479401.83 


479590.65 


171.641 479767.98 


479933.97 


480088.74 


480232.46 


480365.25 


480487.28 


480598.68 


480699.61 


480790.23 


480870.68 


480941.14 


481001.76 


481052.71 


481094.14 


481126.24 


481149.17 


481163.11 


481168.24 


481164.73 


481152.77 


.291 481132.54 


481104.24 


481068.06 


.651 481024.19 


480972.83 


480914.20 


480848.48 


480775.90 


480696.66 


480610.99 


480519.11 


480421.24 


480317.61 


480208.46 


480094.03 


479974.56 


479850.29 


479721.49 


-135.141 479588.40 


479451.30 


-142.641 479310.45 


479166.12 


479018.59 


478868.16 


478715.10 


478559.73 


* V ] 


•M 


• ■MI*] 


552807.89 


549010.29 


545212.69 


541415.09 


537617.49 


533819.89 


530022.29 


526224.69 


522427.09 


518629.49 


514831.89 


511034.29 


507236.69 


503439.09 


499641.49 


495843.89 


492046.29 


488248.69 


484451.09 


480653.49 


476855.89 


473058.29 


469260.69 


465463.09 


461665.49 


457867.89 


454070.29 


450272.69 


446475.09 


442677.49 


438879.89 


435082.29 


431284.69 


427487.09 


423689.49 


419891.89 


416094.29 


412296.69 


408499.09 


404701.49 


400903.89 


397106.29 


393308.69 


389511.09 


385713.49 


381915.89 


378118.29 


374320.69 


370523.09 


366725.49 


362927.89 


359130.29 


355332.69 
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l/AE448data 



378 


sn 


392 


393 


395 


396 


398 


399 


Ena 

PEH 


a 

I 

m 

9 


408 


Em 

eeh 


Em i 


Em 


418 


Em 


370.00 


372.00 


373.00 


374.00 


375.00 


376.00 


377.00 


B 


21810.31 


21960.54 


22112.43 


22266.00 


22421.31 


22578.37 


22737.25 


22897.97 


380 

378.00 

23060.59 

381 

379.00 

23225.14 

382 

380.00 

23391.67 

383 

381.00 

23560.24 

384 

382.00 

23730.89 

385 

383.00 

23903.67 

386 

384.00 

24078.64 

387 

■KEEKEI 

24255.85 

388 

1 386.00 

24435.37 

389 

— 11 >n 

24617.25 


2782997.25 


2804882.54 


2826918.89 


2849107.96 


2871451.47 


2893951.17 


2916608.82 


2939426.28 


2962405.40 


2985548.10 


3008856.33 


3032332.12 


3055977.50 


3079794.60 


3103785.57 




- 158.30 


3176824.16 


389.00 


390.00 


391.00 


392.00 


393.00 


394.00 


395.00 


396.00 


397.00 


398.00 


399.00 


400.00 


401.00 


402.00 


403.00 


404.00 


405.00 


406.00 


407.00 


408.00 


409.00 


410.00 


411.00 


412.00 


413.00 


414.00 


415.00 


416.00 


417.00 


418.00 


419.00 


420.00 


24988.35 


25177.70 


25369.68 


25564.37 


25761.83 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


25962.16 


3226428.09 


3251510.89 


3276784.36 


3302251.16 


3327914.02 


EEEZXEEKB 


3405700.10 




3483586.58 


3509548.74 


3535510.91 


3561473.07 


3587435.23 


3613397.39 


3639359.55 




3691283.87 


3717246.03 


3743208.19 


3769170.36 


3795132.52 


3821094.68 


3847056.84 


3873019.00 


3898981.16 


3924943.32 


3950905.48 


3976867.64 


4002829.81 


4028791.97 


161.08 


161.98 


162.77 


162.65 


162.16 




- 154.04 


147.96 


144.28 


- 130.48 




- 112.33 


- 97.63 


• 120.22 


■ 150.98 


■ 181.74 


■ 212.51 


• 243.27 


■ 274.04 


■ 304.80 


• 335.57 


■ 366.34 


■ 397.11 


■ 427.88 


- 458.65 


- 489.42 


- 520.20 


■ 581.75 


612.53 


643.31 


■ 674.09 


• 704.88 


■ 735.66 


797.24 


828.04 


858.83 


889.63 


478402.33 


478243.23 


478082.74 


477921.18 


477758.89 


477596.20 


477433.46 


477271.02 


477109.26 


476948.54 


476789.23 


476631.74 


476476.45 


476323.79 


476174.16 


476028.01 


475885.76 


475747.87 


475614.82 


475487.07 


475365.11 


475249.45 


475140.62 


475039.13 


474945.54 


474840.70 


474705.10 


474538.74 


474341.61 


474113.72 


473855.07 


473565.64 


473245.46 


472894.50 


472512.78 


472100.28 


471657.02 


471182.98 


470678.17 


470142.58 


469576.22 


468979.08 


468351.16 


467692.46 


467002.98 


466282.71 


465531.65 


464749.80 


463937.16 


463093.73 


462219.50 


351535.09 


347737.49 


343939.89 


340142.29 


336344.69 


332547.09 


328749.49 


324951.89 


321154.29 


317356.69 


313559.09 


309761.49 


305963.89 


302166.29 


298368.69 


294571.09 


290773.49 


286975.89 


283178.29 


279380.69 


275583.09 


271785.49 


267987.89 


264190.29 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 


260392.69 






















































































































































































































































































l / AE 448 data 



XElHini 


85.00 


84.50 


83.00 


234.69 


248.46 


262.45 


291.12 


320.72 


335.89 


382.93 


449.54 


9.5 


9 . 


31 


597.95 


618.07 


78.00 


77.50 


31 


6.5 


6 . 


75.50 


74.50 


724.39 


746.84 


769.70 



80.58 


9.96 


5.6 


Thrust 


0.00 


3635328.12 


3635599.61 


3636068.00 


3636735.59 


3637604.52 


3638676.70 


3639953.82 


3641437.35 


3643128.53 


3645028.31 


3647137.44 


3649456.35 


3651985.22 


3654723.95 


3657672.12 


3660829.04 


3664193.67 


3667764 . 

.69 

3671540 . 

.44 

3675518 

.93 

3679697 

.82 

3684074 

.43 

3688645 

.66 

3693408 

.02 

3698357 

.60 

3703490 

.10 

3708800 

.82 

3714284 

.68 

3719936 

.20 

3725749 

.52 

3731718 

.41 

3737836 

.31 

3744096 

.28 


3750491.08 


3757013.13 


3763654.57 


3770407.26 


3777262.78 


3784212.51 


3791247.58 


3798358.96 


3805537.43 


3812773.65 


3820058.15 


3827381.39 


3834733.78 


3842105.67 


3849487.46 


3856869.55 


3864242.41 


3871596.60 


Pressure 


14.700000 


14.697595 


14.689235 


14.674810 


14.654251 


14.627492 


14.594474 


14.555144 


14.509458 


14.457377 


14.398872 


14.333921 


14.262509 


14.184631 


14.100290 


14.009499 


13.912280 


13.808664 


13.698693 


13.582416 


13.459897 


13.331205 


13.196425 


13.055651 


12.908992 


12.756567 


12.598509 


12.434962 


12.266083 


12.092042 


11.913017 


11.729202 


11.540798 


11.348018 


11.151087 


10.950237 


10.745710 


10.537757 


10.326637 


10.112616 


9.895967 


9.676969 


9.455903 


9.233060 


9.008730 


8.783206 


8.556786 


8.329764 


8.102438 


7.875102 


7.648051 


7.421574 


0.002377 


0.002376 


0.002375 


0.002373 


0.002369 


0.002365 


0.002360 


0.002353 


0.002346 


0.002338 


0.002328 


0.002318 


0.002306 


0.002293 


0.002280 


0.002265 


0.002249 


0.002233 


0.002215 


0.002196 


0.002176 


0.002155 


0.002134 


0.002111 


0.002087 


0.002063 


0.002037 


0.002011 


0.001983 


0.001955 


0.001926 


0.001896 


0.001866 


0.001835 


0.001803 


0.001771 


0.001737 


0.001704 


0.001670 


0.001635 


0.001600 


0.001565 


0.001529 


0.001493 


0.001457 


0.001420 


0.001384 


0.001347 


0.001310 


0.001273 


0.001237 


0.001200 
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l / AE 448 data 


54 

55 

56 I 
5 7 

5 8 

5 9 

6 0 
61 
62 

63 1 

64 

6 5 
6 6 
67 
6 8 

69 r 

7 0 

71 

72 
7 3 
7 4 

75 r 

7 6 
7 7 

78 r 

7 9 
80 | 
81 

82 I 

83 I 

8 4 T 

85 r 

86 r 

87 T 

88 | 

89 [ 

90 r 

92 [ 

93 [ 

94 [ 

95 [ 

96 

97 

98 

99 

100 
101 
102 
103 

’ 10 4 
[ 10 5 
106 


74.00 

73.50 
73.00 

72.50 
72.00 

71.50 
71.00 

70.50 
70.00 

69.50 
69.00 

68.50 

68.00 

67.50 
67.00 

66.50 

66.00 

65.50 
65.00 

64.50 
64.00 

63.50 
63.00 

62.50 
62.00 

61.50 
61.00 

60.50 
60.00 

59.50 
59.00 

58.50 
58.00 

57.50 
57.00 

56.50 
56.00 

55.50 
55.00 

54.50 
54.00 

53.50 
53.00| 

52.50 
52.00 

51.50 
51.0 0j 

50.50 
50.00 

49.50 
49.00 

48.50 
48.00 


_H_ 

816.67 _ 

840.79 _ 

865.34 
890.33 
915.75 
941.62 
967.94 
994.71 
1021.94 
1049.63 ~ 
1077.78 ~ 
1106.40 ~ 
1135.49 
1165.05 
1195.09 

1225.61 

1256.61 
1288.09 
1320.06 
1352.53 
1385.80 
1419.72 
1454.21 

1489.28 
1524.92 
1561.15 

1597.95 
1635.32 

1673.28 ' 
1711.81] 

1750.91 
1790.60 
1830.86 ' 
1871.70~ 

1912.95 
1954.76 
1997.27] 
2040.47| 
2084.35] 

2128.92 
2174.16] 
2220.07 
2266.66 
2313.91 
2361.821 
2410.40| 
2459.631 
2509.52 
2560.06 
2611.25 
2663.10 
2715.59 
2768.73 


I 

73.75 
73.14 ~ 

72.53 
71.92 
71.32 
70.72 
70.12 ~ 

69.53 
68.94 

68.36 
67.77 
67.20 
66.62 
66.06 
65.49 
64.93 

64.37 
63.82 
63.27 
62.73 
62.19 
61.65 
61.12 
60.60 
60.07 

59.56 
59.04 
58.53 
58.03 

57.52 
57.03 

56.53 
56.04] 

55.56 
55.08 
54.60 

54.12 
53.65 
53.18 
52.72 

52.26 
51.80 
51.35 
50.90 
50.45 
50.01 
49.56 

49.13 
48.69 

48.26 
47.83 
47.40 
46.98 


3878922.82] 

3886211.90 
3893454.85 

3900642.90 

3907767.51 
3914820.40 
3921793.57 
3928679.34 
3935470.33 
3942159.56 
3948740.37 

3955206.52 

3961552.16 
3967771.84 
3973860.55 
3979813.71 

3985627.17 
3991297.25 
3996820.67 
4002194.65 
4035115.95 
4045271.09 
4055118.40 

4064655.73 

4073881.73 
4082795.75] 
4091397.87 
4099688.91 
4107670.35 
4115344.34 
4122713.67T 
4129781.75 
4136552.54 
4143030.57 
4134351.il] 
4147546.01 

4160122.37 

4172095.81 

4183481.81 
4194296.37] 
4204555.9l| 
4214277.28 
4223477.63 

4232174.38 
4240385.16 
4248127.74] 
4255419.95 
4262279.66 
4268724.68 
4274772.73 

4280441.39 
4285748.03 
4290709.751 


K _ 

7.195959 _ 

6.971488 _ 

6.748437 _ 

6.527077 

6.307670 

6.090473 

5.875730 

5.663679 

5.454546 

5.248548 

5.045888 

4.846759 

4.651342 

4.459803 

4.272298 

4.088967 

3.909938 

3.735325 

3.565228 

3.399733 

3.238903 

3.082764 

2.931357 

2.784716 

2.642863 

2.505806 

2.373544 

2.246066 

2.123348 

2.005357 

1.892050 

1.783376 

1.679272 

1.579669 

1.484493 

1.393666 

1.307092 

1.224668 

1,146289 

1,071843 

1.001217 

0,934297 

0.870963 

0.811096 

0.754574 

0.701275 

0.651076 

0.603855 

0.559488 

0.517854 

0.478832 

0,442302 

0.408146 


L_ 

0.001163 
0,001127 
0.001091 
0.001055] 
0,001020 
0,000985] 
0.000950 
0,000916 
0.000882] 
0.000849] 
0.000816] 
0.000784] 
0,000752] 
0,000721 
0.000691 
0.00066l] 
0.000632] 
0.0006041 
0.000576] 
0,300550 
0.000524] 
0.000498] 
0.000474] 

0.000450] 
0.000427] 
0.000405] 
0,000384] 
0.000363] 
0.000343] 
0.000324] 
0 000306 
0.000288] 
0 000272] 
0^000255] 
0 !q00240] 
0.000225] 
0 . 000211 ] 
0 loooi9i] 
0.000185] 
0 000173] 
0.000162] 
0 000151] 
0 00014l] 
0.00013l] 

0 000122] 

O.OOOlli] 
0 , 00010 ^ 
0 .ooooii] 

0 .0000901 
0.000084] 
0.0000 77 
0.0000 72 
0 . 0000«>6 
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l / AE 448 data 




4295343.40 


4299665.45 


4303692.05 


4307438.91 


4310921.32 


4314154.13 


4317151.69 


4319927.87 


4322496.04 


4324869.03 


4327059.16 


4329078.22 


4330937.47 


4332647.63 


4334218.91 


4335661.00 


4336983.07 


4338193.79 


4339301.36 


4340313.48 


4341237.42 


4342079.97 


4342847.52 


4343546.04 


4344181.09 


4344757.88 


4345281.22 


4345755.63 


4346185.26 


4346573.98 


4346925.35 


4347242.69 


4347529.02 


4347787.15 


4348019.66 


4348228.91 


4348417.09 


4348586.19 


4348738.07 


4348874.41 


4348996.71 


4349106.36 


4349204.59 


4349292.55 


1739748.39 


1739776.21 


1739800.88 


1739822.77 


1739842.21 


1739859.47 


1739874.80 


1739888.42 


1739900.52 


0.376249 


0.346496 


0.318778 


0.292985 


0.269012 


0.246758 


0.226123 


0.207012 


0.189334 


0.172998 


0.157922 


0.144023 


0.131224 


0.119451 


0.108635 


0.098708 


0.089607 


0.081272 


0.073648 


0.066681 


0.060320 


0.054520 


0.049237 


0.044428 


0.040057 


0.036086 


0.032483 


0.029218 


0.026260 


0.023584 


0.021165 


0.018981 


0.017010 


0.015233 


0.013632 


0.012192 


0.010897 


0.009732 


0.008687 


0.007748 


0.006907 


0.006152 


0.005475 


0.004870 


0.004330 


0.003851 


0.003427 


0.003050 


0.002716 


0.002419 


0.002155 


0.001920 


0.001712 


0.000061 


0.000056 


0.000052 


0.000047 


0.000044 


0.000040 


0.000037 


0.000033 


0.000031 


0.000028 


0.000026 


0.000023 


0.000021 


0.000019 


0.000018 


0.000016 


0.000014 


0.000013 


0.000012 


0.000011 


0.000010 


0.000009 


0.000008 


0.000007 


0.000006 


0.000006 


0.000005 


0.000005 


0.000004 


0.000004 


0.000003 


0.000003 


0.000003 


0.000002 


0.000002 


0.000002 


0.000002 


0.000002 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000001 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 



























































































































































l / AE 448 data 


160 
1611 
162 ~ 
1 631 

I 6 41 
165 
1 6 61 

167 ~ 

168 
169 
17 0 | 
17lf 
17 2 j 
17 31 
17 4 1 ~ 
175) 

17 6 

177 

178 

179 

180 
181 

1 8 2 r 

I 8 3[ 
1 8 41 
185 

II 8 61 

I I 8 7 j 
188 
189 

i 9 pr 

191| 

j 19 2 r 

193 

194 
[1 9 5j 

196] 

I 97 j 

II 9 81 
1991 
200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 

till 

212 


30.001 

30.00 

30.00 

29.50 
29.00 

28.50 
28.00 

27.50 
27.00 

26.50 
26.00 

25.50 
25.00 ~ 

24.50 
24.00 

23.50 
23.00 

22.50 

22.00 

21.50 
21.00 

20.50 

20.00 

19.50 
19.00 

18.50 
18.00 

17.50 
17.00 

16.50 
16.00 

15.50 
15.00 

14.50 
14.00 

13.50 
13.00 

12.50 

12.00 

11.50 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00 
11.00) 
11.00 


H 

6085.45 
6118.12 
6151.04 
6184.23 
6217,70 

6251.45 

6285.47 

6319.77 

6354.35 ~ 
6389.21 ' 

6424.35 ' 

6459.77 ' 

6495.47 

6531.46 ' 
6567 •731 
6604.29 ' 
6641.13 ~ 
6678.26 ~ 
6715.67 ~ 
6753.38 ' 

6791.37 
6829.65 ' 
6868 . 22 ) 
6907.09 ' 
6946.24 ~ 
6985.69 ' 
7025.44 
7065.48 
7105.82 
7146 . 45 ) 

7187.38 
7228.61 
7270 , 14 ' 
7311 . 97 ) 
7354 . 10 ) 
7396 . 53 ) 

7439 . 27 ) 
7482 . 31 ) 
7525 . 66 ' 
7569.31 

7613 . 27 ) 
7657.57 
7702 . 23 ) 
7747 . 25 ) 
7792 . 63 ) 
7838.37 
7884.46 
7930.92 
7977 . 73 ) 
8024 . 89 ] 
8072.41 
8120.28 
8168.51 


28.02 

27.78 

27.53 

27.29 

27.05 

26.81 

26.56 

26.32 
26.07 

25.82 

25.57 

25.32 
25.07 

24.82 

24.57 
24.31 
24.06 
23.80 
23.55 
23.29 
23.03 
22.77 
22.52 
22.25 
21.99 
21.73 
21.47 
21.21 
20.94 
20.68 
20 . 4 l " 
20.15 
19.88 
19.61 
19.35 
19.08 
18.81 
18.54 
18.27 
18.00 
17.73 
17.46 
17 . 19 ' 

16.93 
16.67 

16.41 
16.16 
15.91 
15.66 

15.41 
15.17 

14.93 
14.69 


1739911.29 
1739920.87 
1739929.39 

1739936.97 

1739943.72 

1739949.73 
1739955.09 
1739959.86 
1739964.11 
1739967,90 
1739971.27 

1739974.29 

1739976.97 
1739979.37 
1739981.51 
1739983.42 
1739985.13 

1739986.65 
1739988.01 
1739989.23 
1739990.32 

1739991.29 

1739992.17 
1739992.95 

1739993.65 

1739994.27 
1739994.83 
1739995.34 
1739995.79 
1739996.20 

1739996.56 

1739996.89 

1739997.18 
1739997,45 

1739997.68 

1739997.90 
1739998.09 

1739998.27 
1739998.42 

1739998.56 

1739998.69 
1739998.81 

1739998.91 
1739999.01 
1739999.09 
1739999.17 
1739999.24 
1739999.30 
1739999.36 
1739999.41 
1739999.46 
1739999.50 
1739999.54 


K __ 

0.001527 

0.001362 

0,001215 

0.001085 

0.000969 

0.000865 

0.000773 

0.000691 

0.000618 

0,000552 

0.000494 

0.000443 

0.000396 

0,000355 

0.000318 

0,000285 

0,000256 

0,000230 

0.000206 

0.000185 

0.000167 

0.000150 

0.000135 

0.000121 

0.000109 

0.000099 

0.000089 

0.000080 

0.000072 

0.000065 

0.000059 

0.000054 

0.000049 

0.000044 

0,000040 

0.000036 

0.000033 

0.000030 

0.000027 

0.000025 

0.000023 

0.000021 

0.000019 

0.000017 

0.000016 

0.000014 

0.000013 

0,000012 

0.000011 

0.000010 

0.000009 

0.000009 

0.000008 


L _ 

0.000000 
0.000000 
0.000000 
0,000000 
0.000000 
0.000000 
0,000000 
0.000000 
0,000000 
0.000000 
0.000000 
0,000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0 . 000000 
0.000000 
0.000000 
0.000000 
0.000000 
0,000000 
0.000000 
0.000000 
0.000000 
0,000000 
0.000000 
0,000000 
0.000000 
0,000000 
0.000000 
0.000000 
0,000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0,000000 
0,000000 
0.000000 
0 , 000000 
0.000000 
0,000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
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l / AE 448 data 



G I 

H 

I 1 


_ £ _ ;l 


213 

11.00 

8217.09 

14.45 

1739999.58 

0.000007 

0.000000 

2 14 

11 * 00 

8266.02 

14.22 

1739999.61 

0.000007 

0.000000 

^ X I 

2 15 

11.00 

8315.30 

13.99 

1739999.64 

0.000006 

0.000000 

M A W 

216 

11.00 

8364.93 

13.77 

1739999.67 

0.000006 

0.000000 

4m M V 

2 17 

11.00 

8414.92 

13.54 

1739999.70 

0.000005 

0.000000 

A f 

2 18 

11.00 

8465.25 

13.32 

1739999.72 

0.000005 

0.000000 

2 19 

11.00 

8515.94 

13.10 

1739999.74 

0.000004 

0.000000 

2 2 0 

11.00 

8566.97 

12.88 

1739999.76 

0.000004 

0.000000 

4i V 

2 21 

11.00 

8618.35 

12.67 

1739999.78 

0.000004 

0.000000 

222 

11.00 

8670.08 

12.46 

1739999.79 

0.000004 

0.000000 

223 

11.00 

8722.16 

12.25 

1739999.81 

0.000003 

0 . 000000 

2 2 4 

11.00 

8774.59 

12.04 

1739999.82 

0.000003 

0.000000 

« « ~ 

2 2 5 

11.00 

8827.37 

11.84 

1739999.83 

0.000003 

0.000000 

2 2 6 

11.00 

8880.50 

11.64 

1739999.85 

0.000003 

0.000000 

2 27 

11.00 

8933.98 

11.44 

1739999.86 

0.000002 

0.000000 

228 

11.00 

8987.81 

11.24 

1739999.87 

0.000002 

0.000000 

2 2 9 

11.00 

9041.98 

11.05 

1739999.88 

0.000002 

0.000000 

2 30 

11.00 

9096.51 

10.86 

1739999.88 

0.000002 

0.000000 

2 31 

11.00 

9151.38 

10.67 

1739999.89 

0.000002 

0.000000 

2 32 

11.00 

9206.61 

10.48 

1739999.90 

0.000002 

0.000000 

2 33 

11.00 

9262.19 

10.29 

1739999.90 

0.000002 

0.000000 

2 3 4 

11.00 

9318.12 

10.11 

1739999.91 

0.000002 

0.000000 

2 35 

11.00 

9374.40 

9.93 

1739999.92 

0.000001 

0.000000 

2 3 6 

11.00 

9431.03 

9.75 

1739999.92 

0.000001 

0.000000 

2 37 

11.00 

9488.02 

9.58 

1739999.93 

0.000001 

0.000000 

238 

11.00 

9545.36 

9.40 

1739999.93 

O.OOOOOll 

0.000000 

2 3 9 

11.00 

9603.05 

9.23 

1739999.94 

0.000001 

0.000000 

2 4 0 

11.00 

9661.10 

9.06 

1739999.94 

0.000001 

0.000000 

2 41 

11.00 

9719.51 

8.89 

1739999.94 

0.000001 

0.000000 

2 42 

11.00 

9778.28 

8.73 

1739999.95 

0.000001 

0.000000 

2 43 

11.00 

9837.40 

8.56 

1739999.95 

0.000001 

0.000000 

244 

11.00 

9896.88 

8.40 

1739999.95 

0.000001 

0.000000 

2 45 

11.00 

9956.72 

8.24 

1739999.95 

0.000001 

0.000000 

2 4 6 

11.00 

10016.92 

8.09 

1739999.96 

0.000001 

0.000000 

2 47 

11.00 

10077.49 

7.93 

1739999.96 

0.000001 

0.000000 

2 48 

11.00 

10138.42 

7.78 

1739999.96 

0.000001 

0.000000 

2 4 9 

11.00 

10199.72 

7.63 

1739999.96 

0.000001 

0.000000 

2 50 

11.00 

10261.38 

7.48 

1739999.97 

0.000001 

0.000000 

251 

11.00 

10323.41 

7.33 

1739999.97 

0.000001 

0.000000 

2 52 

11.00 

10385.81 

7.18 

1739999.97 

0.000001 

0.000000 

2 53 

11.00 

10448.58 

7.04 

1739999.97 

0.000001 

0.000000 

2 5 4 

11.00 

10511.73 

6.90 

1739999.97 

0.000000 

0.000000 

2 5 5 

11.00 

10575.24 

6.76 

1739999.97 

0.000000 

0 . 000000 

256 

1 11.00 

10639.14 

6.62 

1739999.97 

0.000000 

0.000000 

257 

11.00 

10703.41 

6.48 

1739999.98 

0.000000 

0.000000 

258 

11 .oc 

10768.06 

6.35 

1739999.98 

0.000000 

0.000000 

25 9 

t 11 .oc 

10833.10 

6.22 

1739999.98 

0.000000 

0.000000 

2 60 

1 11 .oc 

10898.52 

6.09 

1739999.98 

o.oooooc 

0.000000 

2 61 

11 .oc 

10964.32 

5.96 

1739999.98 

o.oooooc 

0.000000 

mm w 4i 

2 62 

] 11 .oc 

11030.52 

5.85 

1739999.98 

0.000000 0.000000 

*2 63 

11 .oc 

11097.1C 

5.7C 

1739999.96 

0.000000 0.000000 

.2 6 4 

11 .oc 

11164.0' 

5.56 

1739999.96 

0.000000 0.000000 

2 65 

11.oc 

11231.44 

5.46 

1739999.96 

0.000000 0.000000 
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>66 


268 


269 


278 


2 

8. 

0_ 

2 _ 

8 

r 

2 _ 

00 

2_ 

2 _ 

CO 

3. 

2 _ 

8, 

£ 

£ 

8 

7 

2 

00 


2 

8 

7 

2 

1 

00 

2 

1 

9 


29 


29 


292 


293 


29 


295 


296 


298 


299 


308 


EEQ 

HU 


312 


gmi 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


318 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


11.00 


H 


11299.21 


11367.38 


11435.95 


11504.92 


11574.30 


11644.10 


11714.30 


11784.92 


11855.96 


11927.43 


11999.31 


12071.63 


12144.38 


12217.56 


12291.18 


12365.24 


12439.75 


12514;71 


12590.12 


12665.99 


12742.32 


12819.11 


12896.38 


12974.11 


13052.33 


13131.02 


13210.21 


13289.88 


13370.05 


13450.72 


13531.90 


13613.59 


13695.79 


13778.52 


13861.77 


13945.56 


14029.88 


14114.75 


14200.17 


14286.14 


14372.68 


14459.78 


14547.46 


14635.72 


14724.57 


14814.02 


14904.06 


14994.72 
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Figure 3.0.2: Plot of X-Position vs. Time 
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Figure 3.0.4: Plot of Y-Position vs. Time 
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Figure 3.0.5: Plot of Weight vs. Time 
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Figure 3.0.6: Plot of Inclination vs. Time 
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Figure 3.0.7: Plot of Total Velocity vs. Time 
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Figure 3.0.9: Plot of Thrust vs. Time 
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Figure 3.0.10: Plot of Thrust/Weight Ratio vs. Time 
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Figure 3.0.11: Plot of Y-Position vs. X-Position 
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4.0 Orbital Analysis 

The NASA ground rules for the Shuttle II specified two orbits into which the orbiter must be 
capable of launching: a 270 nautical mile orbit due east of Kennedy Space Center, and a 150 
nautical mile polar orbit To verify the feasibility of the proposed ASTS configuration, it 
was necessary to determine the types and magnitudes of the orbital maneuvers that would 
be necessary for orbit inserton. As a starting point, the output data for the final position and 
velocity from the launch simulation program was used. This data places the orbiter at an 
initial altitude of 78.1 nautical miles with a velocity of 25,962 feet per second. 


4.1 The Insertion Orbit 

The initial position and velocity of the orbiter place the orbiter originally at the perigee of 
an elliptical orbit. The semi major axis of this orbit, ai, was calculated using the Vis-Viva 

equation: . ^ . 

v ‘- (n ‘ ^ 

where ri is the sum of the altitude of the orbit and the radius of the earth (20925672.57 feet.). 

The eccentricity of this first orbit, and the altitude and velocity at the apogee of this orbit 
were calculated using simple geometric relationships. The pertinent data on the first orbit 

is as follows: 



ei = .02473 
ai = 21,943,244 feet 
riap = 22,485,871 feet 
Viap = 24,709 feet/second 


4.2 The 270 Nautical Mile Orbit 

In order to provide a minimum energy transfer between the insertion orbit and the 270 nmi 
orbit, an elliptical transfer orbit was used (see figure 4.2). 

The 270 nmi orbit has a radius of &2 = 22,567,272 ft. The semi major axis of the transfer 
ellipse, at, has a value of (ai +&2)f2 - 22,255,258 ft. 

The velocities at periapsis and apoapsis of the transfer ellipse, as well as the velocity for 
the circular orbit were calculated using the Vis*Viva equation again. 









The following values were calculated for the 
insertion into a 270 nmi orbit: 

Vlap 

Vti 

dVi 


velocities and the velocity changes needed for 

24,709 ft/sec 
24,890 ft/sec 

181 ft/sec 


V t 2 24,799 ft/sec 

V 270 24,975 ft/sec 

dV2 176 ft/sec 

% 

The total impulsice change in velocity for this orbital maneuver is 351 ft/sec. 

4.3 The 150 Nautical Mile Orbit 

The path of the elliptical orbit transfer orbit used for this maneuver is shown in figure 4.3. 

The velocity and radius of the 150 nmi orbit are found to be 

r 150 21837672ft 

Viso ’ 25389 ft/sec 

Following the sane procedure as before in order to obtain the velocity changes necessary 
yields the following information: 

Vlap 
Vtl 

dVi 

V t2 
Vi50 

dV 2 


24709 ft/sec 
24678 ft/sec 

-31 ft/sec 

25419 ft/sec 
25389 ft/sec 

-30 ft/sec 


The total impulsive change in velocity for this orbital transfer is 
velocity changes would have to be made in a direction opposite to 
the shuttle. 


-61 ft/sec. These two 
the direction of travel of 
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Figure 4.3.1:150 Nautical Mile Polar Orbital Transfer 
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5.0 Stability and Control Analysis 

Stability and Control Analysis for the Advanced Space Transportation System was 
performed using geometric Shuttle data and empirical and graphical methods given by 
Etkin, Roskam and Smetana. The geometric parameters used to perform the analysis 
included physical data for the wing, the body, and the vertical tail. 

In order to determine stability and control derivatives, the geometric data for the Shuttle II 
must be used to determine aerodynamic parameters such as aspect ratio, wing taper ratio, 
and the mean aerodynamic chord. 

'A' is the aspect ratio and is defined to be equal to t^/S; b is the wing span and S is the wing 
planform area. 'S' can be easily determined using the following formula: S=(c r +ct)(b/2), 
where c r is the wing root chord and ct is the wing tip chord. The ratio of the wing tip chord to 
the wing root chord is in itself and important physical parameter known simply as the 
taper ratio, X. The value of the mean aerodynamic chord (the wing chord which passes 
through the centroid of either half-wing) can be found by integrating the equation for the 
square of the wing chord: c 2 = (c r - [(c r -ct)/(b/2)]y) 2 over the span of half a wing and 
multiplying by 2/S. All of these aerodynamic parameters have been calculated and are 
found in the Executive Sum mary in Tables 3.2.1 and 3 j.l. 

5.1 Determination of Cia 

CLa is an extremely important stability derivative describing the change in lift of an 
aircraft due to variation in angle of attack. 

Roskam gives the following formula for determining CLa : 

Clo = Kwh^Low 

where K w b =1 - .25(d/b) 2 + 0.25(d/b) 

and CLaw = 2jtA/(2+[(A 2 B 2 /K 2 )(l+tan 2 Ac/2/B 2 )+4] 1/2 } 

In these equations: d = body diameter at root chord = 30 ft. 

b = wing span = 127.7 ft. 

A = wing aspect ratio = 2.8 

P = compressibility parameter = (1 - M 2 ) 172 (for subsonic flight) 

K =s ratio of actual wing section left curve slope to 2 n 
Ac/ 2 = sweep of the wing half-chord line = 34° 

Since these equations are useful only for Mach numbers less than the critical Mach 
number, K is assumed to be approximately 1. Hence, the following values for Cl<x are 

obtained: 


M range G 2 

0.0 - 0.3 1.0 

0.3 - 0.6 0.775 

0 .6-Mcr 0.415 


f fc. (/rad) 
2.97 
3.09 
3.33 


It is apparent that these values alone for the lift-curve slope (although correct) cannot be 
used to represent CLa for the ASTS. For this reason, methods of estimating CLa at higher 
Mach numbers and varying Reynolds numbers (from Dynamics of Flight, by Etkin) were 
employed. 
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Etkin presents a graphical method by which the lift-curve slope of an aircraft may be 
obtained as a function of its geometrical parameters. Before incorporating Etkin s method, 
values for t/c and r (defined in Figure 5.1.1) were assumed to be 0.06 and 5 , respectively. 

Using Figure 5.1.1 the following section lift curve slopes were obtained: 


Ba 


lajlgi/fcadl 

IQ? 

0.87 

5.74 

10 7 

0.93 

6.14 

10 8 

0.95 

6.27 


It follows that these values for the actual lift curve slope are good approximations for 
transition both at the leading edge and at c/2. 

These values of (ai)o were then used in conjunction with Figure 5.1.2 to estimate the 
overall aircraft lift curve slope. Also, Etkin defines K as [(ai)om/2*]P and Ap as the 
compressibility sweep parameter equal to tan' 1 (tanAc/Vp). 


Table 5.1.1. Clo of ASTS for M < Mcr 


BCiVK </°) I C La (/rad) 



Note that these values of Clu are slightly higher than those estimated using the method 
outlined by Roskam. 


Etkin presents a system for predicting Cl« for Mach numbers greater than the critical 
Mach number using the graphical method shown in Figure 5.1.3 which is a function o 
M, X, and Ac/2. Results using this figure are listed in Table 5.1.2. 
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Table 5.1.2. CLa of ASTS for M > Mo- 


Mrange 

AtanAc /2 

A(M z -lW* 

CLa (/rad) 

1.0 -1.3 

1.89 

1.64 

4.14 

1.3 -1.6 

1.89 

2.97 

3.42 

1.6 - 1.9 

L89 

4.04 

2.72 

1.9 - 2.2 

L89 

5.02 

2.21 

2.2 - 2.5 

L89] 

5.97 

1.82 

2.5 - 2.8 

1.89 

r 6.88 

1.62 


Prediction of CLa for Mach numbers greater than approximately 2.8 was not possible due to 
the lack of a proper estimation technique. 

5.2 Determination of the Neutral Point, hn 


The overall aircraft neutral point, hn, is defined as the point where the pitching moment is 
invariant with angle of attack. It is further defined geometrically by Figure 5.2.1. 

Roskam gives the following relationship for determining the neutral point . 

Tin = hnw + Ahnn (for tail*less configurations) 


where hnw = ki[(XacVc r )-k2] 

and Ahnfi = - (dM/da)/(qScCLaw) 

where dM/da = (q/36.5) S Wf^CXilde/daAXi 


Note that Roskam defines both the neutral point, hn, and the aerodynamic center, Xac/c, as 
the same quantity. 


Using the relationships for Ahn B , it was found that the body had little or no effect on the 
overall aircraft neutral point Therefore, determination of the neutral point is reduced to 

the following form: 


hn = ki[(Xac'/c r )-k 2 ] 

Xac'/cr ki and ko are found using Figures 5.2.2, 5.2.3, and 5.2.4, respectively. Using 
these figures together with known geometric data and the relationship for hn, the following 

results were obtained. 


Table 5J2.1. Neutral Point of the ASTS for M < Mo- 
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These values for hn show that as the Mach is increased, the neutral point begins to move 
gradually aft. This finding agrees with most aerodynamic data on similar aircr 


It is important to determine the characteristics of a flight vehicle in as many ofits 
regimes as possible. Therefore, methods given by Etkin were used to find values of h 
Mach numbers greater than the critical Mach number. 


Etkin gives a graphical method shown in Figure 5.2.5 for determining hn as a function of 
A, M, 3W and Ad2- Using these variables and Figure 5.2.5, the following data was obtained. 


Table 5A1 Neutral Point of the ASTS for 0 < M < 2& 


M range 

Atan Ac /2 

Ad-M 2 ) 1 ^ 

hn 

0.0 - 0.3 

1.9 

2.8 

0.28 

0.3 - 0.6 

1.9 

2.5 

0.29 

0.6 - Mcr 

1.9 

1.8" 

0.31 

1.0 - 1.3 

1.9 

1.6 

0.44 

L3 -1.6 

1.9 

3.0 

0.'14 

1.6 -1.9 

1.9 

4.0 

0.48 

1.9 - 2.2 

1.9 

5.61 

0.48 

2.2 - 2.5 

1.9 

6.0 

0.49 

2.5 - 2.8 

1.9 

6.9 

0.49 


i 



Again, these numbers show that (as expected) the neutral point starts out relatively close to 
the quarter-chord for low subsonic Mach numbers and moves aft approaching the ha 
chord at higher supersonic Mach numbers. 


5.3 Determination of Cmo 

C ma is an extremely important stability derivative which dictates whether or not a 
configuration has positive or negative static stability. Cma is determined through the 

following equation: 


Cma = ^La (h ■ hn) 


where h represents the non-dimensional location of the aircraft center of gravity. For a 
configuration to be statically stable, Cma must be negative; this requires that the quantity 
(h - hn) must be negative. Xcg was determined by considering the location and weight of 
the components of the ASTS. The distance between the center of gravity and the leading 
edge of the mean aerodynamic chord (non-dimensionalized by dividing by the mean 
aerodynamic chord) is defined as h. Upon completion of these calculations, it is found that 

h = 0.136. 


Hence, Cma = Clxx (0.136 - hn) 



5.4 Clq of Vertical Tail 
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The lift curve slope of the vertical tail is determined by the geometric characteristics of the 
vertical tail and is influenced by the change in flow characteristics due to thepresenceof 
the wing and body combination. Figure 5.4.1 shows the relationship between (ap)B/ai for a 
vertical tail on a body of circular cross section. Using a value of 1 for D/h, (aF)B/ai is 
approximately 2.9. 


Values ofcSfor the vertical tail are determined using the same procedure as before for the 
wing. 'WhWthis value is obtained it is multiplied by the factor (aF)B/ai to correct for its 
location mthe flow field relative the the wing and body. It is assumed that the vertical tail 
has approximately the same values of t/c and x as the wing so that both the wing and 
vertical tail have the same values for (ai)o. Using the geometric parameters for the 
vertical tail and Figures 5.1.2 and 5.1.3 the following is found. 



Table 5.4.1. Corrected Clo values for the ASTS Vertical Tail 



Determination of Cut corrected for the vertical tail beyond a Mach number of 
approximately 1.6 was not possible. 


5.5 Determination of C| P 

Ci„ is the change in rolling moment due to a variation in rolling velocity. Ci p is often 
referred to as the rolling moment damping derivative. Etkin provides a graphical method 
by which Ci p can be calculated based on the following parameters: p, K r A,^»a nd Ae (which 
is defined as tan* 1 [(l/pltanAc/^). Etkin’s method is presented here in Figure 5.5.1. Table 
5.5.1 contains the values of C] p found using Etkin s method. 
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Table 5.5.1. Q p values ofASTS for M < Mcr 



Cip does not vary greatly with Mach number, and is always negative. 


5.6 Determination of Cip 

Cib represents the change in rolling moment due to a change in sideslip angle, p. Qp is 
primarily governed by the wing dihedral angle, the body and the vertical tail. For this 
reason, Cip is often referred to as the effective dihedral derivative. Eltkin presents a 
graphical method of determining Qp based on the effects of the wing, the body and the 
vertical tail. This method utilizes Figure 5.6.1. Using this figure, it was found that 
Cip w = *.32 Cl. 


The effect of the body on Cip is found using the equation: 


(ACip)B = 1-2 (A) 1/2 Zw (h+wl/b 2 

where Zw = vertical distance of c/4 point below fuselage centerline = 11.0 ft. 
h s average fuselage height at wing root = 28.0 ft. 
w = average fuselage width at wing root * 30.0 ft. 


Using this information: (AC]p)B = 0.08. 

The effect of the vertical tail on Cip is found using the equation. 


(ACip)p = -ap (1 - do/dpXVp/V)2 Sp Zf/(S b) 

where Zp = vertical distance between vertical tail AC and aircraft CG = 41.0 ft. 
Sp = Vertical Tail planform area = 603 ft 2 . 

Using this equation under the assumption that da/dp is negligible and Vf/V is 
approximately equal to 1: (ACip)p = -0.033 ap 


Combining components: C]p= -0.32 Cl + 0.08 - 0.033ap 

This equation holds for .all Mach numbers below the critical Mach number. 
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5.7 Determination of Ci r 

Clr represents the change in rolling moment due to yawing velocity, or simply rolling 
moment due to yaw. Etkin divides Ci r into two components: Ci r due to the wing and Q r due 

to the vertical tail. Using Figure 5.7.1, (Ci r ) B = 0-23 C L- 

The contribution of the vertical tail on Ci r is shown through the following formula: 

(Qr )F = S F Z F t* +(lf/b)]/(S b) 

where If * distance between vertical tail AC and the aircraft CG = 48.0 ft. 

Hence, (Q r )f ■ 0.079 aF 

Combining components: Ci r = 0.23 Cl + 0.079 aF 

This equation holds for all Mach numbers below the critical Mach number. 

5.8 Determination of C n p 

CnB is defined to represent the change in yawing moment due to sideslip angle, |J. C n p is a 
very important stability derivative since it alone determines whether or not an aircraft has 
positive weathercock stability; for positive weathercock stability, C n p must be positive. 

Figure 5.8.1 shows the variation of C n p/C L 2 as a ftmction of wing aspect ratio, A, and wing 
quarter-chord sweep angle, A. 

Using Figure 5.8.1, C n p = 0.07 Cl 2 

5.9 Determination of C n p 

C np is a cross derivative defined as the change in yawing moment due to rolling velocity. 
Etkin presents the following equation for determining Cn p : 

Cnp = [(ACnp)l/ClJ Cl + [(AC np > 2 /(CD 0 )a] (Cdo)<* 
where (AC np )i/CL and (AC np )2/(CD0)a are found from Figure 5.9.1 
From Figure 5.9.1, (AC np )i/CL = -0.17 and (AC np )2/(CDo)<x = 19 , so that C np may be written 


as 


Cnp = -0.17 Cl + 19 (Cdo)cl 


5.10 Determination of Cnr 

Cnr is defined to be the change in yawing moment due to yawing velocity and is often 
referred to as the damping in yaw derivative. Etkin presents the following equation for 
determining C n r: 

Cnr = [(AC nr )i/CL 2 ]CL 2 + KACmO^CdoICdo 
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where [(AC^i/Cl 2 ] and [(ACm-WCdo! are found from Figure 5.10.1 

Figure 5.10.1 reveals that [(ACnr)l/CL 2 ] = -0.02 and [(ACnr)2/CD0] = -0.48, so that C nr may be 
written as 

Cm- = -0.02 C L 2 * 0-48 QDO 

5.11 Determination of CMq 

C mq is defined as the variation in pitching moment due to pitch velocity or simply as the 
damping in pitch. Smetana gives the following relationship for C mq : 

C mq = - 2 x’ lx*l Clo/c 2 

where x’ = distance from CG to wing quarter-chord (positive for forward CG) = -6.97 ft. 

After substituting the proper geometric and stability parameters, C mq can be written as 
Cmq — 0.0372 Cj^a 

5.12 Determination of Ciq 

CLq is defined as the change in lift coefficient due to pitch velocity. Smetana gives the 
following relationship for CLq • 


CLq = 2 x Clo/c 

Upon substitution of known geometric parameters, CLq Is written as 
CLq = -0.2727 Clcc 

5.13 Determination of Cd<x 

Cdo is defined as the change in drag coefficient due to a change in the angle of attack. 
Smetana gives the following relationship for CDa: 

Cdc = (dCuo/da) + [2 Cl/(*A.)] CLa 

Since (dCD 0 /da) is usually small, Cdo may be written simply as [2 Cl/(xA)] Clo- 


Therefore, Cd« may be expressed as 
Cdo = 0.2274 Cl Clo- 

5.14 Determination of C y p 

Cyft is defined as the change in side force,y, due to a change in sideslip angle, |J. Smetana 
breaks Cyp up into 3 components: that of the wing, body, and vertical tail. 

(Cyp)w = Cl 2 (6 tan A sin A)/[A(A + 4 cosA)] 

Upon substitution, (Cyp) w = 0.2232 Cl 2 
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(C y p)B = -ki (CLa)B (Sb/S) 

where Sb = body reference area = 3733 ft* 

ki = interference factor determined from Figure 5.14.1 


Figure 5.14.1 gives ki as a function of Zw/(d/2) where Zw is the distance from 0 y 
centerline to quarter-chord point of exposed wing root chord (positive for quarter-chord 
point below body centerline) and d is the max body height at wing/body intersection. 


Zw = H O ft 
d = 25.0 ft 


From Figure 5.14.1, ki = 1.44 

Smetana gives the following values of (Cl<x)b : 

(Clo) circular body = 0.0525 rad* 1 
(Clo) rectangular body = 0.1253 rad* 1 

Since the ASTS configuration is a combination of these two shapes, let (Clo)b = 0.08 rad 1 

Now that all values needed are known, (Cyp)fl = -0.074 

% 

(Cyp)F = -k (Clo)F d + da/dp) Sp/S 
k is found from Figure 5.14.2 and is a function of [bV(2ri)l. 


by = vertical tail span = 29.0 ft. 

ri = fu se lage diameter in tail area = 30 ft 


From Figure 5.14.2, k = 0.75 

If it is assumed that do/dB is negligible then (Cyp)F may be expressed as 

(Cyp)p s -0.078 aF 

Combining all components; Cyp * 0.2232 Cl^ - 0.074 - 0.078 aF 

5.15 Determination of Cyp 

Cjrp represents the change in side force due to rolling velocity. Smetana expresses C yp as 
the sum of the contributions due to the wing and the vertical tail. 

(Cyp ) w = Cl ([(A + cosA) tanMA + 4 cosA)] + 1/A) 

Substitution into (Cyp )w yields (Cyp )w = 0.91 Cl 


(Cyp )f = -2 aF Sf Zp/(S b) 

All of these parameters have been defined previously so that (Cyp )F = -0.067 aF 
Combining these two; Cyp = 0.91 Cl - 0.067 aF 
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5.16 Determination of Cyr 

Cyr is defined as the change in side force due to yawing velocity. Smetana expresses Cyr 
asthe sum of the contributions due to the wing and the vertical tail. 


(Cyr)w = 0.143 Cl • 0.05 

(Cyr)P = -2 If (Cyr)F/b = -0.752 (Cyr)F 



Once the proper stability derivatives have been determined, it is possible to use them to 
investigate the dynamic stability of the configuration. Presented here are the longitudinal 
and lateral stability matrices corresponding to the non-dimensional, controls fixed, 
linearized equations of motion for a rigid aircraft. 


Longitudinal Eouations: 

Dx = Ax where x = (u a q 0) T 

A(l,l) = (Cxu + 2 Clo tan0o)/2n 

A(l,2) = CxoAi 

A(l,3) = 0.0 

A(l,4) * -Clo/2p 

A(2,l) * (Czu - 2 CloV(2ji - Czo) 

A(2,2) = Czo/(2|i - Czo) 

A(2,3) = (2p. + CZqV(2M. - Czo) 

M2,4) = -Clo tan©o/(2|i - Czo) 

A(3,l) = (C M u + Cmo A(2,l))/iB 
M3£) = (CMa + C mo A(2,2)ViB 
A(3,3) = (CMq + Cma A(2,3)yiB 
A(3,4) = Cma A(2,4)/iB 
A(4,l) = 0.0 
A(4,2) = 0.0 
A(4,3) = 1.0 
A(4,4) = 0.0 

where p. = m/(p S c/2) 

ifi = Iyy/[p S (c/2) 3 ] 

Cxa = Clo - CDo 
Cza = -CLa * Cdo 
Czu = - (M 2 /l-M 2 ) c lo 
CZq ~ -CLq 
CDm = Cdu/M 

Cxu = *2 (Cdo + Clo tan©o) * M Com 
Cza = 0.0 
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T.flt.ftral Equations: 

Dy = By where y = (B p r O)™ 

B(l,l) = C y p/2p 
B(l,2) — Cyp/2p. 

B(l,3) = -(2p - CyrV2p 
B(l,4) = CL0/2H 
B(2,l) = GC Qp + iE C n pyA 
B(2,2) = (ic Ci p + iE C np )/A 
B(2,3) = (ic Clr + iE Cn^A 
B(2,4) = 0.0 

B(3,l) = ftA C n p + iE QpYA 

B(3,2) = (iA C np + iE Cl p )/A 

B(3,3) = (iA Cnr + iE ClrYA 

B(3,4) = 0.0 

B(4,l) = 0.0 

B(4,2) = 1.0 

B(4,3) = tan0 O 

B(4,4) = 0.0 


where p = m/(p S h/2) 

A = iAiC-iE 2 
iA = Ixx/[p S (h/2)^] 
ic = Izz/[p S (h/2)^] 
ig = Ixz/[p S (b/2)^] 

Once the matrices A and B are found, their Eigen Values can be determined. These Eigen 
values reflect whether or not the system is stable. For a system to asymptotically stable, all 
of its Eigen values must have negative real parts. 

The matrices A and B were calculated using a computer program written in FORTRAN. 
The program and results of the investigation of the dynamic stability for landing 
(assuming 0 O = 10° and U 0 = 300 ft/sec and L = W) are listed on pages 66 and 71. 



5.18 Determination of Ci&j 


Ci 5 a describes the variation in rolling moment with change in aileron deflection; often, 
this stability derivative is referred to as the aileron power. Smetana gives the following 
equation for Ci5a- 

Ab 

Ci6a = [2CLaX/(Sb)]\ cydy 

where c is the equation for the wing chord as a ftinction of wing span, y. 


c = c r - [(c r -cty(h/2)] y = 73.125 - 0.863 y 

and a and b represent the spanwise distance for the body centerline to the inner and outer 
aileron points, respectively. 



a = 27.0 ft. b = 58.5 ft 

x is a correction factor which is a function of the aileron chord to wing chord ratio, x is 
determined from Figure 5.18.1. Since the ratio of aileron chord to wing chord is not 
constant over the span of the aileron, the average value was used to determine x. 

(c a /c w ) av g = 0.33 fcs0.57 

Upon integration and substitution of the proper values, 

Ci5o = 0.071 Clo. 

5.19 Determination of Cisr 

Ci5r is the variation in rolling moment coefficient with change in rudder deflection. 
Smetana gives the following formula for estimating C]5r • 

ClSr = ap t Sp Zp/(S b) 

Here, t is a function of the rudderarea to vertical tail area ratio. 

Sr = 157.5,fjfT| Sp=603ft2 

Hence, Sj/Sp = 0.281 so that from Figure 5.19.1, t = 0.48. 

Substitution into the equation for Cigr results in 
Cisr s 0.0158 ap. 

5.20 Determination of C n sr 

Cn8r is the variation in yawing moment coefficient with a change in rudder deflection and 
is often referred to as the rudder power. Smetana gives the following relationship for C n 6r 
C n 5r = -aF t Sf lf/(S b) 

As with Cisr, the value oft is determined from Figure 5.19.1 and has a value of 0.48 so that 
Cngr may be expressed as 

CnSr * *0.0187 ap 
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AE448 CG Locations/I 



A ] 

B 

c 

D 

1 

Component 

X-Location 

Weiaht 

Y-Moment 

2 

wina 

116.60 

28,300 

3 , 299,780 

3 

Vertical Tail 

140.00 

3,200 

448,000 

4 

Forward Body 

28.90 

27,000 

780,300 

5 

Mid Body 

80.00 

23,600 

1 . 888,000 

6 

Aft Body 

128.00 

13,200 

1 , 689,600 

7 

Main Enaines 

137.00 

22,800 

3 , 123,600 

8 

OMS 

138.00 

3,700 

510,600 

9 

RCS 

138.00 

1,600 

220,800 

EQ 

Control Systems 

48.00 

4,000 

192,000 

WWW 

Avionics 

22.00 

r 3,100 

68,200 

WWW 

Personnel 

35.00 

1,600 

56,000 

WWW 

Stores 

4 2 . 00 

600 

25,200 

m 

Encrine Accessories 

125.00 

I ” 6,200 

775,000 

WWW 

Hvdroaen Tanks 

64.00 

24,000 

1 , 536,000 

m 

Oxyqen Tank 

121.00 

11,000 

1 , 331,000 

EH 

TPS 

96.30 

r 6,100 

587,430 


Totals 


180,000 

16 , 531,510 

EQ 

Centroids 

91.84 
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*%% 


CGNTROI S FIXED LONGITUDINAL.. AMD I..A1 f-.lv hL 
DYNAMIC STABILITY AND CONTROL PROGRAM 
FOR THE SECOND GENERATION SPACE 
TRANSPORTATION SYSTEM * 

WRITTEN BY* PHILIP BENEFIELD % DEL JOHNSON 


C 


DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUBLE 

DOUELE 

DOUBLE 

DOUBLE 


PRECISION 

PRECISION 

PRECISION 

PRECISION 

PRECISION 

PRECISION 

PRECISION 

PRECISION 

PRECISION 


Q* S V w* us * CBAR* IJO * RHO * THET AO 

B*A*C *E *CLO*CDO 7 CMU 1 CMA 

CMAD*CMQ*CLA*CLO *CDA * CDU 

CLB y CLP y CLR *CNB* CNP y CNR 

CYB •! CYP y CYR y LLONG y MUL.ONG » IBL.ONG 

CX A y CZA y SOUND y M y CZIJ y CZG y CDM 

CXU y CZADy ALONG(4 y 4)y LLAT 

IALATyICLATyIELATy MULAT 

DEL y ALAT < 4 y 4) 


0PEN(UNIT=9y FILE='STAB1IN* DAT'y ST ATUS= ' OLD') 

OPEN< UNIT-10 y FILE='STAB1OUT.DAT'ySTATUS®'UNKNOWN 

READ ! 9 y & ) G y S y U y US y CBAR y UO y RHO y THET AO 

READ(9 y #) By Ay Cy Ey CL0 y CD0 

READ(9 y *) CMU y CMA y CMAD y CMQ 

READ<9y*> CLAyCLQyCDAyCDU 

•READ! 9 y*) CLB y CLP y CLR 

READ(9 y %) CNB » CNP y CNR y CYB y CYP * CYR 



C 

300 

20 

% 


FORMAT <IX y A40 y' = 'yF16.6) 

FORMAT <////////y10Xy'GEOMETRICALyINERTIAL 

' characteristics: '»//> 


X EQUILIBRIUM 


y 


c 



WRITE<10y20) 
URITE<10 y 300) 
WRITE < 10 y 300 ) 
WRITE(10 y 300) 
WRITE!10y 300) 
WRITE!10f300) 
WRITE!10 y 300) 

WRITE!10y300) 


'GRAVITATIONAL ACCELERATION !FT/SEC~2) ' • 0 
'WING AREA !FT~2 )'yS 
'WEIGHT !LB)'yU 
'WING SPAN !FT)'yWS 

'MEAN AERODYNAMIC CHORD !FT)'yCBAR 

'AIRSPEED !FT/SEC)'yUO 

'AIR DENSITY ! 5LUGS/FT~3)'*RHO 


WRITE(10y300) 

WRITE!10y300) 
WRITE!10y300) 
WRITE!10y300) 

WRITE!10y300) 
WRITE!10 y 300) 
WRITE!10y300) 


'INITIAL THETA !RAD)'yTHETAO 


FORMAT!////* 


WRITE!10y30) 

WRITE!10 y 300) 
WRITE!10*300) 
WRITE!10*300) 
WRITE!10*300) 
WRITE!10y300) 
WRITE!10*300) 
WRITE!10y300) 
WRITE!10y300) 'CDU 


'IVY ! 

SLUG-FT-2)'»B 

'IXX ! 

SLUG-FT~2) ' yA 

'IZZ (SLUG-FT 2) yC 

' IXZ !SLUG-FT”2)'*E 

'CLO 

(-) ' * CLO 

'CDO 

(-) ' y CDO 

X , ' NOh 

J-DTMENSIGNAL 

' CMU 

!-) ' y CMU 

'CMA 

!-) ' y CMA 

'CMAD 

!-) ' y CMAD 

'CMQ 

!-) ' * CMQ 

' CL A 

!-) ' y CL A 

' CLQ 

!-) ' y CLQ 

'CDA 

(-) ' , CDA 


_ STABILITY DFRIOATT 


m rr c ♦ ' , 


PAGE IS 

Of poor 


nriii 


o n n 


40 


FOR M A Tv// / / '> 10 X y ' L A T E R A L-DIRECTIO N A L 


WRITE(10 •/ 40) 
WRITE(10 y300) 
WRITE(10 y 300) 
WRITE(10*300) 
WRITE(10 y 300) 
WRITE<10?300) 
WRITE <10 y 300) 
WRITE(10 y 300 > 
WRITE(10 y 300) 
WRITE(10 y 300) 


' CLB </RAD)'yCLB 
'CLP </RAD)'yCLP 
'CLR (/RAD)'y CLR 
'CMP (/RAD)'y CNB 
'CNP (/RAD)'»CNP 
'CNR </RAD)'yCNR 
'CYB (/RAD)'y CYB 
'CYP </RAD)'y CYP 
'CYR (/RAD)'yCYR 


DERIVATIVES X 


/) 


%%% CALCULATIONS FOR LONGITUDINAL MOTION 


LL0NG=CBAR/2♦0 

MULONG=(W/G)/(RHO*S*LLONG) 

IBL0NG=B/(RH0*S*<LL0NG**3)) 

CXA=CLO~CDA 

CZA=~CLA“CDO 

S0UND=1116.4 

M=UO/SOUND 

CZU=-( (M#M) / ( 1 »0-M*M) ) fcCLO 

CZQ=-CLQ 

rriM-rnu/M 

CXU=~ < 2♦0#(CDO + CLO#TAN CTHETAO)))-<M*CDM) 

CZAD=0.0 

ALONG(1y1) = (CXU+2.0*CLO*TAN < THETAO))/< 2.0*MUL0NG) 
ALONG<1y 2)=CXA/(2♦O&MULONG) 

ALONG<ly3)=0»0 

ALONG<1y4)=-CLO/(2 *0*MUL0NG) 

AlONG(2 y1) = (CZU-2.0*CL0)/(2♦0*MUL0NG-CZAD) 

ALONG< 2y 2)=CZA/(2♦0*MUL0NG-CZAD) 

AL0NG<2y3) = (2» 0#MUL0NG+CZQ)/< 2♦OKMULONG CZAD) 

ALONG(2 y 4) = (-CLO*TAN(THETAO))/C 2.0*MUL0NG-CZAD) 

ALONG (3 y 1) = ( CMUPCMAD*ALONG ( 2 1 1) )/IBLONG 

ALONG< 3y 2) = (CMA+CMAD#ALONG< 2y 2))/IBLONG 

ALONG v 3 y 3) = <CMQ+CMAD*ALONG(2 y 3))/IBLONG 

ALONG < 3 y 4) =CMAD*ALONG ( 2y 4) /'IBLONG 

ALONG(4y1)*0.0 

ALONG< 4y2)=0.0 

ALONG(4 v 3)=1»0 

ALONG <4 y 4 )= 0.0 


C 


v \V 

/pv 




CALCULATIONS FOR LATERAL MOTION 


M JLAT = < W/G) / ( RHO*S*LLAT) 
IALAT = A/< RH0*S#LLAT**3) 
ICLAT=C/<RHC#S*LLAT#*3) 
IELAT-E/ <RH0&S*LLAT**3) 
DEL=(IALAT#ICLAT)-IELAT*#2 


C 



c 


10 

c 


15 

C 

100 


101 


200 

201 

u 


ALAT (1,1 )=CYB/<2.0#MULAT) 

AL AT(1 9 2>-CYP/(2 * OKMULAT) 

ALAT (1,3 ) - - ( 2 * 0#MULAT-CYR ) / < 2. O'-KiiULAT ) 
ALAT(1,4)-CLO/<2♦0*MULAT) 

ALAT(2y1)«<ICLAT#CU3+IELAT*CNB)/DEL 
ALAT(2,2)=(ICLAT*CLP+IELAT#CNP)/DEL 
ALAT( 2t 3)*(IC L A T % C L E + IELATKCNR)/DEL 
ALAT(2y 4)=0♦ 0 


ALAT(3?1)=<IALAT#CNB+IELAT*CLB)/DEL 
Al AT(3 , 2) = (IALAT*CNP+IELAT*CLP)/DEL 
ALAT(3y3)*(IALAT#CNR+IELAT*CLR)/DEL 


ALAT (3,4)=0*0 
ALAT(4, 1)=0»0 
ALAT(4 ,2)=1 ♦ 0 
ALAT < 4 ,3)=TAN(THETAO) 
“ SLWt ? 4T4T=0‘ rtf- 


WRITE(10 y 100 ) 

DO 10 J = 1y 4 

WRITE <10 y 200) (ALONG(J y K)y K=ly4) 

IF < J♦NE♦4) WRITE <10y201) 

CONTINUE 


WRITE(10 y101) 

DO 15 J=l,4 

WRITE(10 y 200) (ALAT(J y K)y K=ly4> 

IF(J♦NE.4) WRITE (10*201) 
CONTINUE 


* 

* 

* 

* 

* 

* 


FORMAT(//////*25X y'THE STABILITY MATRIX CORRESPONDING TO',/, 
?OX,'TWE NONDIMENSIONALy CONTROLS FIXED, LINEARIZED,'»/, 
25Xy'LONGITUDINAL EQUATIONS OF MOTION FOR'y/, 


35X y ' A RIGID AIRPLANE ' y//) _ P , 

r ORMAT < //// y 25X y ' THE STABILITY MATRIX CORRESPONDING TO 
20Xy'THE NONDIMENSIONALy CONTROLS FIXED, LINEARIZED,' 
23X,'LATERAL EQUATIONS OF MOTION FCR'y/y 
35X,'A RIGID AIRPLANE',//) 


-ORMAT(5X,'5 ' y4<F12.3,6X)*' !') 


y / y 
/ 

f / t 


FORMAT 


rr' / 
w /v 


y / 


74X 


) 


CLOSE(9) 
CLOSE(10) 

ST CP 
END 
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GEOMETRICAL * INERTIAL & EQUILIBRIUM CHARALI ERlRTxC..' <• 


GRAVITATIONAL ACCELERATION (FT/SEC 2) 

WING AREA (FT~2) - 
WEIGHT (LB) = 
WING SPAN (FT) = 
MEAN AERODYNAMIC CHORD (FT) = 

AIRSPEED (FT/SEC) = 
AIR DENSITY (5LUGS/FT~3) * 
INITIAL THETA (RAD) * 
IYY <SLUG-FT~2) = 
IXX (SLUG-FT~2) = 
IZZ (SLUG-FT"2) = 
IXZ (SLUG-FT"2) = 

CLO (-) = 

CDO (-) = 


32♦174000 

5818.000000 
180000.000000 
127.692000 
51.120000 
300.000000 
0.002377 
0.174500 
12490393.000000 
1686164,000000 
13399266.000000 
795369.000000 
0.418000 
0.050000 


NON-DIMENSIONAL STABILITY 

CMU 

CMA 

CMAD 

CMQ 

CLA 

CLQ 

CDA 

CDU 


derivatives: 

(_) = 0.000000 

(_) = -0.403200 

<_) = 0.000000 

(_) = 0.104200 

( _> = 2.800000 

<_) = -0.763500 

( _) s 0.266000 

( _ 1 = 0.000000 


lateral-directional 


derivatives: 


CLB 

(/RAD) 

— 

- 0,338000 

CLP 

(/RAD) 

sr 

- 0.225000 

CLR 

(/RAD) 

= 

0.775500 

CNB 

(/RAD) 


0.012000 

CNP 

(/RAD) 


- 0.070 000 

CNR 

(/RAD) 

=: 

- 0.030000 

CYB 

(/RAD) 


- 0,706000 

GYP 

(/RAD) 

— 

- 0.195800 

CYR 

(/RAD) 

— 

0,51420 0 



-0*00315899 
-0.02743637 
0.00000000 
0.00000000 


-0.05570903 

-0.74066202 

-0.04074188 

0.00000000 


r i i r~ 

n CL 


THE STABILITY MATRIX CORRESPONDING 10 
N 0 N DIM E N S10 N A L ? CONTROLS FIXE D •> LINE A R1 L h. D * 
LONGITUDINAL EQUATIONS OF MOTION FOR 

A RIGID AIRPLANE 


0.00480166 

-0.09003117 

-0.00745438 

0.00000000 


0.00000000 
1.02411888 
0.00192645 
1.00000000 


'\ /v I ^ /\ /; nr 
kj * 0 1 o ‘t ■ 


-0.00232787 


0.00000000 

0.00000000 


THE STABILITY MATRIX CORRESPONDING TO ^ 

THE NONDIMENSIONAL i CONTROLS FIXED» LINEAR! /.ED? 
LATERAL EQUATIONS OF MOTION FOR 
A RIGID AIRPLANE 


-0.01545018 

-0.50320992 

-0.04867205 

1.00000000 


-0.95942552 

1.69903754 

0.09279544 

0.17629303 


0.03298353 

0.00000000 

0.00000000 

0.00000000 


70 





• ¥V« 


lies* 'ja-a* Cw ( -i.^T a T * £cs--- 

uaS' ■ *- j - s » w • • — - * ■» - •• ^ * - • • 

1 44 4 S >« 4 toS 4* •* 4'' I ^ , ■-•*•4— * .*.» • • i k J 


LINEAR SYSTEM ANALYSIS RESULTS 
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Figure 5.1.1. Lift curve slopes for two-dimensional incompressible flow. 
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Figure 5.1.2.Lift-curve slop, for swept and tapered wings at speeds below critical Mach 
number. 
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Figure 5.1.3. Lift-curve slope for swept and tapered wings at supersonic speeds. 
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Figure 5.2.2. Aerodynamic Center Location of Lifting Surfaces. 
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Figure 5.2.4. Aerodynamic Center Transformation Constant K 2 
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Figure 5.2.5. Chordwise position of the mean aerodynamic center of swept and tapered 
wines at high speeds expressed as a fraction of the mean aerodynamic 


chord. 







Figure 5.4.1. Lift-curve slope for single fin and rudder on a 


body of circular cross section 








Figure 5.5.1. C] p for straight-tapered wings 
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Figure 5.6.1. Qp for straight-tapered wings with no dihedral. 
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Figure 5.7.1. Charts 


for estimating Ci r for subsonic incompressible flow. 


A, degrees A, degrees 







, . . c .2 with aspect ratio and sweep for subsonic incompressible 

Figure 5.8.1. Variation of C n p/CL wtn asp 
flow 
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Figure 5.10.1. Charts 


for estimating (AC nr )l/CL 2 a« d (AC nr)2/CD0 


a £ 








Figure 5.14.1. 


Values for wing-fuselage interference factor. 
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Figure 5.14.2. Values for k as a function of the ratio of vertical tail span to fuselage 

disinter in the tail region. 
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Figure 5.18.1.Values for x as a function of aileron chord to wing chord ratio. 



Figure 


5.19.1. Values for t hs a function of rudder area 


to vertical tail area ratio. 




